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Abstract 

Thla   pap)er,  on  the  diffraction  of  waves  by  smooth  convex  obstacles,  con- 
sists of  two  parts.  In  Fart  I  a  new  geometrical  theory  of  wave  propagation 
is  explained  and  applied  to  the  problem.  In  Part  II  the  diffraction  problem 
for  various  obstacles  is  formiJ.ated  as  a  boiindary  value  problem  and  solved 
by  the  usual  separation  of  variables  method,  Hhe   solutions  are  then  asymp- 
totically expanded  for  large  ka  =  2na/\  (a  »  obstacle  dimension,  X  "  wove- 
lene;th)  and  compared  with  the  results  of  Part  I.  In  all  cases  the  geometric 
solution  is  found  to  agree  with  the  asymptotic  expansion  of  the  exact 
solution. 

The  geometrical  theory  of  Part  I  is  developed  for  solutions  of  the  time 
reduced  Maxwell's  equations  as  well  as  the  reduced  wave  equation.  The 
diffracting  objects  are  assumed  to  be  of  arbitrary  convex  shape  and  of  any 
tjrpe  (e.g.,  acoustically  hard  or  soft,  perfect  conductors,  dielectrics). 
In  addition  to  developing  and  applying  the  geometrical  theory  to  variously 
shaped  bodies  in  Part  I,  a  certain  feature  of  the  self-consistency  of  the 
theory  is  demonstrated.  Soma  calculations  are   also  carried  out  to  indicate 
the  numerical  usefulness  of  the  geometric  theory. 

The  exact  solution  of  the  diffraction  problem  associated  with  the  reduced 
wave  equation  or  the  time  reduced  Maxwell's  equations  for  the  bodies  considered 
in  Part  I  are  found  and  asymptotioa^ly  expanded  in  Part  II.   The  results  of 
this  section  are  an  extension  of  previous  work  in  the  types  of  boxmdary  condi- 
tions and  problems  considered.  The  problem  of  asymtitotioally  expanding  the 
electromagnetic  field  diffracted  by  an  imperfectly  conducting  circular  cylinder 
is  treated  as  well  as  other  vector  and  scalar  diffraction  problems  associated 
with  spheres  and  cylinders. 
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Part  I.  Application  of  the  Geometrical  Theory  of  Diffraction 
1,   Introduction 

This  paper,  on  the  diffraction  of  a  wave  by  a  Errooth  convex  opaque  object  of 
any  shape,  consists  of  two  parts.  In  Part  I  the  "  peometrical  theory  of  diffrac- 
tion" L  J  J  a  new  theory  of  wave  propagation,  is  explained  and  applied  to  the  prob- 
lem. In  Part  II  the  diffraction  problem  for  various  objects  of  special  shape  is 
formulated  as  a  boundary  value  problem  and  solved  by  the  usual  separation  of 
variables  method.  The  solutions  are  then  expanded  asymntotically  for  large  ka  ■ 
— Y"  (a  -  obstacle  dimension,  \  •>  wavelength)  and  compared  with  the  results  of 
Part  I.  In  each  case  the  leading  terms  in  the  expansion  of  the  solution  coincide 
with  the  simple  geometrical  results,  thus  verifying  the  geometrical  theory. 

The  •'  geometrical  theory  of  diffraction"  is  an  extension  of  geometrical  optics 
which  is  also  based  upon  the  postulate  that  fields  propagate  along  rays.  However, 
it  introduces  new  kinds  of  rays,  called  diffracted  rays,  in  addition  to  the  usual 
optical  ones.  In  the  present  case  these  diffracted  rays  lie  in  part  on  the  surface 
of  the  diffracting  object,  so  they  may  also  be  called  surface  rays.  Furthermore, 
the  new  theory  assigns  a  field  value,  which  includes  a  phase,  an  amplitude, and, 
in  the  electrcm.agnetic  case,  a  direction  or  polarization  to  each  point  on  a  ray. 
The  total  field  at  a  point  is  postulated  to  be  the  sum  of  the  fields  on  all  rays 
which  pass  through  the  point.  By  using  this  theory  we  will  obtain  an  explicit  ex- 
pression for  the  field  produced  at  any  point  when  a  wave  hits  a  smooth  convex  opaque 
object. 

The  theory  is  developed  for  both  scalar  and  vector  (e,g,,  electromagnetic) 
fields  and  for  objects  of  any  type  (e.g,,  acoustically  hard  or  soft,  perfect  elec- 
trical conductors,  dielectrics).  The  results  depend  upon  the  nature  of  the  object 
in  an  essential  way.  To  simplify  the  analysis  we  will  ignore  rays  which  go  through 
the  object,  thereby  restricting  our  results  to  opaque  objects. 

An  Important  advantage  bf  the  present  method  is  that  it  does  not  depend  upon 
separation  of  variables  or  any  similar  procedure.  Therefore, the  shapes  of  the 
objects  to  which  it  can  be  applied  are  not  restricted.  But  for  simplicity  we  have 
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limited  our  considerations  to  smooth  convex  objects.  These  include  cylinders  of 
convex  cross  section,  which  have  previously  been  treated  by  this  same  method^-'. 

From  its  sjirdlarity  to  geometrical  optics,  the  present  method  can  be  expected 
to  yield  good  results  when  the  wavelength  X  is  small  compared  to  the  obstacle  di- 
mensions. As  has  been  pointed  out  above,  this  is  verified  in  Part  II  in  certain 
cases  where  other  descriptions  of  the  field  are  available  for  comparison.  This 
method  has  also  been  applied  to  diffraction  by  an  aperture  of  any  shape  in  a  thin 
screen  •- -I ,  to  diffraction  in  waveguides  L^-i  and  to  other  problems  L-'J,  in  those 
cases  as  well  as  in  the  present  one,  it  has  been  found  that  the  results  are  useful 
even  for  wavelengths  as  large  as  the  relevant  dimensions  of  the  scatterer. 

Most  of  our  theo3*y  is  independent  of  the  nature  of  the  field  under  consideration. 
Therefore, the  consequences  of  this  part  of  the  theory  are  the  same  for  all  fields  to 
which  the  theory  applies.  For  example, the  procedure  for  locating  the  rays  is  the 
same  for  alD  fields.  On  the  other  hand,  some  aspects  of  the  theory  depend  speci- 
fically on  the  kind  of  field  and  the  type  of  obstacle.  Consequently  these  aspects 
must  be  considered  separately  in  each  case.  In  particular, the  decay  rates  asso- 
ciated with  the  field  on  a  siu-face  ray,  the  diffraction  coefficients,  and  the 
polarization  effectsare  in  this  category.  Therefore, it  will  be  convenient  first 
to  present  the  theory  for  a  scalar  field  and  then  to  modify  the  analysis  where 
necessary,  in  order  to  treat  vector  fieldSg 

The  problem  which  we  consider  is  that  of  finding  the  time  harmonic  field  which 
results  when  a  given  incident  time  harmonic  field  impinges  upon  a  sniocth  object  in 
a  homogeneous  medium.  The  incident  field  will  be  characterized  in  terms  of  the  in- 
cident rays  and  the  field  on  these  rays, 

2o   Ray  tracing 

According  to  geometrical  optics,  the  only  rays  which  occur  in  the  present 
problem  are  the  incident  rays,  the  reflected  and  refracted  rays  which  they  produce 
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on  hitting  the  objent,  and  the  resulting  multiply- reflected  and  refracted  rays. 
All  of  these  rays  ar^  straight  Ijnes  and  their  directions  are  detenrdned  by  the 
laws  cf  reflection  and  refraction.  We  will  ignore  the  refracted  rays  since  we 
assume  the  body  to  be  opaque.  Then,  since  the  body  is  convex,  no  multiply- reflected 
rays  can  occur. 

In  the  geometrical  theory  of  diffraction, additional  rays,  called  diffracted 
rays,  are  introduced.  In  the  present  problem  these  rays  are  produced  by  incident 
rays  which  are  tangent  to  the  surface  of  the  body  [see  Figure  i\,     Fach  tangent 
ray  "  splits"  at  the  point  of  tangency,,  One  part  continues  along  the  path  of  the 
incident  ray.  Another  part  travels  along  the  surface  of  the  body  as  a  surface  ray. 
This  surface  ray  ia  a  geodesic  or  shortest  path  on  th^  body  surface.  At  each  point 
it  splits  again,  one  part  continuing  along  the  surface  ray  and  the  other  part  leaving 
the  surfa3e  along  the  tangent  to  the  surface  rayo  Thus  a  single  grazing  incident 
ray  gives  rise  to  infinitely  many  diffracted  rays,  one  of  which  is  shed  at  each 
point  of  the  surface  ray. 

Diffracted  rays  can  be  characterized  by  an  extended  form  of  Fermat's  Principle, 
In  the  present  case  the  relevant  extension  is  this:  The  diffracted  rays  connecting 
two  points  Q  and  P  in  the  exterior  of  the  body  are  those  curves  joining  Q  and  P  which 
have  stationary  (optical)  length  among  all  curves  from  Q  to  P  having  an  arc  on  the 
body  surface.  From  this  principle  it  follows  that  a  diffracted  ray  consists  of  a 
straight  line  segment  from  Q  to  the  body  surface,  a  geodesic  arc  along  the  surface 
and  a  straight  line  segment  from  the  bod^'  to  Po  It  also  follows  that  the  straight 
line  segments  are  tangent  to  the  arc  These  are  just  the  properties  of  diffracted 
rays  which  were  described  above. 

Diffracted  rays  penetrate  the  usual  shadow  regions  of  geometrical  optics  and 
account  for  the  non-zero  fields  there.  They  also  modiiy  the  field  in  geometrically 
lit  regions. 


-  3a 


>^ 

H- 

^ 

*1 

■» 

fVJ 

•  • 

n 

-d 

•H 

•~i 

o 

D- 

H- 

b-i. 

<•              — 

D 

3 

n 

C* 

O 

H- 

►^ 

T» 

O 

H- 

OJ 

• 

D 

;— • 

1 

M 

T3 

Q) 

0> 

>!• 

3 

H-" 

H- 

O 

1-^ 

» 

■I 

O 

C+ 

Q- 

"-> 

»• 

t— 

(D 

H- 

9 

fi 

O 

•^ 

Q) 

►-^ 

<: 

< 

tu 

» 

o 

c;- 

^-> 

rit 

-I 

1 

T 

o 

<i 

ill 

3 

05 

c+ 

C+ 

c-t- 

0) 

c 

v 

1 

"-i 

lU 

a 

O 

-( 

c 

1i 

o 

Oo 

H^ 

rr 

T5 

* 

d- 

TJ 

T 

3- 

0; 

o 

05 

r-J 

►t 

M 

s: 

a 

.■5 

jj 

»j 

o 

o 

M 

« 

-^ 

C 

5 

+ 

■-» 

3 

M 

c* 

D> 

» 

c* 

3 

a 

V 

O- 

►1 

r»- 

o 

•o 

3* 

c 

M 

(• 

TO 
3- 

+ 

0> 

CO 

>-^ 

o 

01 

>• 

•1 

3 

^ 

o 

1 

01 

3 

c+ 

ft> 

I-* 

rr 

(->• 

3 

T) 

n 

?n 

TO 

J 

3- 

M 

w 

r+ 

.T) 

c+ 

c 

O 

1 

M 

Q- 

t-ti 

!-"• 

H- 

-t) 

a> 

3 

M) 

• 

o 

(B 

"-^ 

9 

►1 

C* 

a> 

•-^ 

05 

o 

1 

3 

c+ 

o 

CD 

(t 

3 

3 

a- 

£) 

«+ 

►< 

M 

0) 

cf 

«< 

<+ 

O 

O 

>-^ 

rt- 

►1 

•n 

3- 

o 

M 

-» 

3 

^ 

rn 

lO 

3 

^ 

c+ 

CL 

0) 

O 

0} 

o 

TJ 

<» 

• 

W 

ti- 

05 

M 

c+ 

Hi 

0) 

3- 

H- 

£) 

a 

a<i 

M 

3- 

^ 

>i 

r+ 

o> 

a> 

Vi 

O. 

H- 

:^ 

o 

►^ 

(0 

o 

»-b 

o 

3 

1 

a. 

0) 

0) 

« 

H- 

O 

Cfl 

CO 

(:•■ 

H- 

n- 

(1) 

O 

0) 

CL 

0) 

o 

1 

•i 

H> 

0) 

O 

«< 

B> 

-  h  - 

3.   Field  calculation 

In  the  geometrical  theory  of  diffraction,  just  as  in  ordinary  geometrical 
optics,  a  field  value  is  associated  with  each  ray.  This  field  value  if  composed 
of  a  phase  function  0(s)  and  an  amplitude  A(s),  both  of  which  may  vary  with  distance 
s  along  the  ray.  For  vector  fields  there  is  also  an  orientation  or  polarization 
which  we  will  treat  later  by  permitting  A(s)  to  be  a  vector.  The  phase  <j(   is  deter- 
mined by  the  optical  law  d^/ds  ■  1  from  which  it  foil  owe  that  0(s)  ■  0+8, 
Here  0^  is  the  phase  at  the  point  Q  from  which  s  is  measured;  moreover,  s  increases 
in  the  direction  of  propagation. 

The  amplitude  A(s)  is  determined  by  the  optical  form  of  the  principle  of  con- 
servation of  energy.  This  principle  states  that  in  time  harmonic  fields  the  energy 
flux  is  the  same  thr«u,<?;h  every  cross  section  of  a  tube  of  rays.  For  a  narrow  tube, 
having  cross-sectional  area  do"  at  a  point  P,  this  flux  is  proportional  to  A  dc. 
Here  A  is  the  amplitude  of  the  field  at  the  point  P  of  a  ray  contained  in  the  tube. 

If  A  and  d<r  are  the  amplitude  and  cross -sectional  area  at  another  point  Q  of  the 

2      2  1/2 

same  ray,  then  we  have  A  dcT  »  A  dot.   Thus  A  •  A  (dcT/dcr)  '  ,  To  evaluate  the 

area  ratio  we  observe  that  dC   is  Just  the  area  cut  out  by  the  tube  on  that  wave- 
front  ^  •  constant  which  passes  through  P.  Similarly^ d<r  is  the  area  cut  out  on 
the  wavefrcnt  0  -  0^  through  Q.  If  s  denotes  the  distance  from  P  to  Q^a  simple 
computation  shows  that  \see   Figure  f]  since  the  rays  are  straight  lines. 


(1) 


^^  .     "1^2 


Here  p^  and  p^  are  the  principal  radii  of  curvature  of  the  wavefront  through  Q, 

Upon  combining  the  above  results  for  ^  and  A,ve  find  that  the  field  u(P) 
associated  with  the  straight  ray  from  Q  to  P  is 


(2) 


u(P)  -  A, 


P1P2 


(p^+8)(p2+s) 


1/2   ik(0^+s) 
e 
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Here  k  ■  w/c  is  thfi  propagation  constant,  co  is  the  angular  frequency  of  the  field j 
and  c  is  its  propagation  speed.  We  see  that  the  field  at  P  is  determined  in  terms 
of  the  field  u(Q)  •  A^e     at  Q.  The  incident  field  Uj^(P)  is  of  the  form  (2) 
since  the  incident  rays  are  straight  lines. 

To  deteririne  the  field  on  a  reflected  ray  we  assume  that  at  the  point  of  re- 
flection this  field  is  proportional  to  the  incident  field.  The  coefficient  of  pro- 
portionality is  the  reflection  coefficient  R  which  depends  upon  the  nature  of  the 
field,  the  properties  of  the  reflector  at  the  point  of  reflection,  the  angle  of 
incidence  and,  in  some  cases,  the  wavelength  X  ■  2tiA«  P"or  vector  fields  R  is  a 
matrix.  Because  R  depends  only  upon  properties  at  the  point  of  reflection,  it  can 
be  determined  from  the  analysis  of  the  simpler  problem  of  reflection  of  a  plane 
wave  from  an  infinite  plane  interface .  Since  R  is  known  for  many  types  of  fields 
and  reflectors, we  will  not  consider  it  further.  The  calculation  of  the  radii  of 
curvature  of  the  reflected  wavefront  from  the  knowledge  of  the  incident  wavefront 

r6i 

and  the  reflector  shape  has  also   been  performed  before"-  -■.  Therefore, we  may  assume 

that  the  reflected  field  u  (?)  can  be  calculated  completely.  It  will  also  be  of 

the  form  (2)  since  the  reflected  rays  are  also  straight  lines. 

We  denote  by  u  (?)  the  geometrical  optics  field  at  Fo  This  is  the  sum  of  the 

fields  -  each  of  the  form  (2)  -  on  the  various  incident  and  reflected  rays  through 

Po  In  the  shadow  u  -  0  since  there  are  no  incident  or  reflected  rays  thereo 

g 

Let  us  now  apply  the  preceding  considerations  to  a  diffracted  ray  such  as  that 

shown  in  Figure  1,  Such  a  ray  consists  of  a  straight  incident  ray  through  Q  tangent 

to  the  surface  at  Q. ,  a  geodesic  arc  on  the  surface  from  Q,  to  ?_,  and  a  straight 

diffracted  ray  from  P^  to  P.  The  field  on  the  incident  ray  up  to  the  point  Q-  is 

^  ik9f^(Q^)    ^ 

gi\en  by  eouation  (2)  abo\'B.  ¥e  now  assume  that  u,(Q-)  ■  A,(Q,  )e       ,  the  field 

^  ^    ^  ^  ik^^CQl) 

on  the  diffracted  (or  surface)  ray  at  Q,,  is  proportional  to  u.(Q,)  =•  A.(Q,)e       j 
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the  field  on  the  Incident  ray  at  Q^    .     We  also  assume  that   the  phases  are  equal 
so  5fj(Q-))  "  0j(Q-i).     Thus  the  proportionality  becomes 

(3)  A^(Q^)      -     r(Q^)A^(Q^)    . 

We  call  the  coefficient  of  proportionality  D(Q, )  the  diffraction  coefficient.  We 
assume  that  it  depends  upon  the  nature  of  the  field,  the  properties  of  the  object 
at  Q,  and  the  propagation  constant  k.  In  the  case  of  vector  fields  D(Q, )  is  a 
matrix. 

Before  attempting  to  determine  the  diffraction  coefficient  D(Q, ),let  us  con- 
sider how  the  diffracted  field  varies  along  the  diffracted  (or  surface)  ray.  If 
t  denotes  distance  along  the  diffracted  ray  from  Q.  then  the  previous  considerations 
yield, for  the  phase  0.(t)jthfl  result  ^A^)   '  ^AQ-,)   +  to  To  determine  the  ampli- 
tude Aj(t)  we  apply  the  energy  conservation  principle  to  a  narrow  strip  of  diffract,ed 
rays  on  the  surface  of  the  object  (see  Figure  3).  Let  us  denote  the  width  of  the 
strip  at  distance  t  from  0^  by  dott)o  Then  the  energy  flux  through  a  cross  section 
of  the  strip  is  proportional  to  A,(t)  drrtt).  At  the  slightly  greater  distance  t  +  dt 
the  flux  will  be  smaller  because  diffracted  rays  are  shed  in  the  inter^/al  from  t  to 
t  ♦  dt. 

The  shedding  of  these  rays  corresponds  to  radiation  of  energy  from  the  strip 
of  surface  rays.  We  assume  that  the  radiated  energy  is  proportional  to  A.(t)  and 
to  the  area  d<5'(t)dt  from  which  radiation  occurs.  Thus  the  energy  conservation 
principle  becomes 

(h)      A^(t-fdt)dr(t+dt)  -  A^(t)dOtt)  -  -2a(t)  A^(t)d(5'(t)dto 


The  case  in  which  the  field,  or  some  component  of  it,  is  required  by  a  boundary 


condition  to  be  zero  on  the  surface  will  be  considered  separately  later.  The 
results  in  that  case  are 
different  interpretation. 


results  in  that  case  are  the  same  as  those  described  here,  but  u,  then  has  a 


Lgure  ?:     A  baad  of  diffracted  rays  on  the  surface  of  s  diffracting 

rbiect.     The  geodesic  distance  frcr  the  point  Q^  is  t,  while 
~«r    is  the  width  of  the  band. 


P(r,fl) 


Q(/?,o) 


The  two  faad-llee  cf  rays  erianating  frcpi  a  line  sotirce  locs'^sd 
at  S  and  diffracted  arcts-;-*   a  circular  cylinder  tc  a  point  ?. 
The  Tipper  faacLly  is  incident  on  the  cylinder  at  Q^    and  lea^^s 
Ire    cylinder  taneentially  at  F,,  while  the  lower  fanily  is  in- 
cnent  on  the  cylinder  at  ^^   and  leaves  the  cylinder  tangential" 
ly  at  ?.. 
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Here  2a'!t)   is  a  prcporticroality  factor  relating  the  ciange  of  flux  *. :    ..-.^  c-^rsy 
radiated.     The  rhvsical  nature  of   :r.=  raiiaticn  rrccers  sugres^s  -.r.a:.   z  : ,    iezer.: 
only  upon  local  properties  cf  -.-.-    r.rfa^e   ar.l  -.he  field.     Vs  irill  ccnsider   ■_'.= 
determination  of  a(t)   la^er.        :-    -e   liviie    '_,    ct  z:  is.t  let   t;  vend  *.c   :-?r:, 
obtaining 

«)  It  [*^^^]  •  -"[4  ^'^]  • 

Upen  integrating  (5)  we  find 

/"zir  \i/2 


(6) 


Here  dc /dc  is   the  rati:   of  -.he  strip  width  at  t  «   I   tc  that  at  t,   cr  rather  tbe 
liisit  cf  this  rati:   as  the  strip  width  taids  tc   zer:  .     Tre  pcint  t  ■  C  is  the 
point  of  ta-a-enr:^  :,    so  k^{2)   z  ^jjC^)  is  siven  by     •      =c.-.    :r    -57  :e    called 

r  ;..    .      I_-    ve  insert  this  value  into   (6^,  replace  d<r  by  d«?(3- .^,   3.--   ir:lide 
the  phase,   ve  cbtairj  f:r  the  diff 


(7)  .,r 


L-      J  1-  J   .„  J 


It  is  t:   te  ncted  that  u.it^:    is  rrcrcrtiCEal  -.  :   --:.r   :.r  -  tent  field  at  ;,  , 

C  "       "  -. 

Tc  obtain  the  diffract-ed  field  at   a  pcint   ?  off  the  s--rfa:s   ■e  first   aptly    J,, 

This  yields  the  diffracted  field  at  the  pcint  :f  tangency  P^  frcn  which  a  tangent 

ray  travels  tc  ?.     Then  we  saJce  -ase  of  (2),  which  applies  tc  a  straight  ray.   tc 

describe  the  variation  of  the  field  frcr  ?.    tc  ?.     Hcvevsr, we  carj-.ct   aprly  (2) 

directly  since  one  cf  the  radii  of  the   diffracted  w^avefrcnt  -   say  p„  -     is  rerc 

at  P,  0     This  is   s:   oecaise  the  rurface   cf  the  cb-^ct  is  a   caustic  cf  the  faailv 
1  "■  • 

of  diffracted  rays,     Tnerefcre  we  proceed  as  fcllrvs; 

Let  Q  dencte   a  pcint  en  the  diffracted  ray  frcr  ?,   tc  ?,      Ihen   ,2,    expresses 
the  field  at  ?  in  terM  of  that  at  ;.     Xcw  let   ;  t«nd  tc  ?,  .     Th^n  o.  t^nds  to 
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zero,  p^  tends  to  some  finite  liMt,  0    tends  to  0A?-,)  and  s  tends  to  the 
distance  from  P^^  to  P  while  u(P)  remains  constant.     It  follows  from   (2)  that 
A    must  become  infinite  in  such  a  way  that  A„  pZ'     has  a  finite  limit.     If  we 


o  ^2 


denote  this  limit  by  A^(P^)  then  (2)  becomes 
(8) 


The  expression  (8)  for  u,(P)  becomes  infinite  at  the  surface  of  the  object, 
at  which  s  «  0  -  a  consequence  of  the  fact  that  this  surface  is  a  caaptic  of  the 
diffracted  rays.  This  is  an  instance  of  the  failure  of  the  present  theory  to  yield 
the  correct  result  for  the  field  at  ar  near  a  caustic.  As  a  consequence  the  result 
must  be  appropriately  modified  in  such  places,  and  we  will  show  how  later.  At  pre- 
sent, however,  this  fact  prevents  us  from  identifying;  the  field  given  by  (8)  with 
that  given  by  (7),  when  both  are  evaluated  at  P, »  If  such  an  identification  were 
possible  we  would  be  able  to  determine  the  amplitude  factor  A  (P  )  in  (8).  In- 
stead we  will  determine  this  factor  by  observing  that  it  must  be  proportional  to  the 


amplitude  in  (7)  evaluated  at  P..  We  write  this  proportionality  as 

1 


(9) 


A^(Pj)  -  D(P^)D(Qj)Aj^(Q^) 


'dctQ^) 
doTP-L) 


r   t 


exp 


i(t)dr 


We  call  the  proportionality  factor  P(P,  )  a  diffraction  coefficient.  It  depends 
upon  the  nature  of  the  field  and  the  object  in  the  neighborhood  of  P^,  We  further 
assume  that  D(P, )  is  the  same  finction  of  these  variables  as  is  DCQ^)*  This  assump- 
tion is  based  upon  the  r eciprocity  principle  -  that  a  source  at  Q  produces  the  same 

field  at  P  as  a  source  at  P  produces  at  Q,  With  oiu"  assumption  this  principle  is 
satisfied,  but  not  otherwise.  Upon  inserting  (9)  into  (8)  we  obtain  for  the 


diffracted  field  at  P  the  result 
(10) 


u^(P)  -  A^(Q^)'T)(P3^)D(Q3^) 


fBtf(Q^) 


dffTPp" 


[sTp^]  «xp  Jik^^(Qi).t*s]-j  a(t)dt 
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In  the  derivation  of  this  result  we  excluded  cases  in  which  the  field  or  some 

component  of  it  was  required,  by  a  boundary  condition,  to  vanish  on  the  surface.  This 

was  necessary  because  our  analysis  dealt  with  the  function  A,(t),  the  diffracted 

field  amplitude  on  the  surface.  However, as  we  have  seen,  the  surface  is  a  caustic 

of  the  diffracted  field.  From  this  it  follows  that  the  field  is  much  stronger  in 

a  surface  layer  some  few  wai-elengths  thick  than  it  is  at  points  a  large  number  of 

wavelengths  from  the  surface.  Therefore, the  preceding  argument  will  still  apply 
and  include  the  previously  excluded  cases  if  A,(t)  is  any  measure  of  the  field 

amplitude  in  this  layer.  The  field  within  the  layer  will  be  characterized  by  a 
definite  profile  along  a  surface  normal,  i.e.,  a  definite  variation  with  distance 
from  the  surface.  The  amplitude  of  the  field  at  any  point  other  than  a  zero  of 
the  profile  car  serve  as  the  measure  of  the  field  amplitude  in  this  layer. 

In  general  the  field  within  the  caustic  layer  is  composed  of  a  number  of  different 
modes,  each  with  its  own  profile.  Each  mode  is  characterized  by  its  own  amclitude, 

A,  and  therefore  its  own  diffraction  coefficient  D  (0,)  and  proportionality  factor 
dm  m  X        - 

for  radiation  a  (t).  Consequently  the  diffracted  field  (lO),  which  represents  the 
m 


contribution  of  only  one  mode,  must  be  replaced  by 


■d<5-(Q,) 


1 
?  r  P,   -I? 


(11).,(F)  .  A.(Q^)exp  Ak[0,(Q^).t*s]l  |^^J   [^^] 

^ 

Equation  (11)  yields  the  diffracted  field  associated  with  a  single  diffracted 
ray,    Jf  more  than  one  ray  passes  through  P,  then  u,(P)   is  a  sum  of  terras  of  the 
above  form,  one  term  for  each  diffracted  ray  through  P.     The  total  field  at  P  is 
then 

(12)  u(P)     -     Ug(p)     *     u^(P). 
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The  above  result  enables  us  to  calculate  the  field  at  any  point  which  does 

not  lie  on  a  caustic  or  shadow  boundary.  This  calculation  involves  the  incident 

field,  various  geometrical  quantities,  the  reflection  coefficient  R,  the  diffraction 

coefficients  D  (Q, )  and  the  decay  exponents  a  .  Our  result  is  not  yet  complete  be- 
n  X  m 

cause  we  have  not  shown  how  to  determine  D  and  a  •  Ihis  is  done  in  the  next  section* 

IB     n 

We  will  later  show  how  to  calculate  the  field  on  or  near  the  surface,  which  is  a 
caustic,  and  at  other  caustics. 


h»       Diffraction  by  a  circular  cylinder 

From  (3),  which  defines  the  diffraction  coefficient  D  (Q, ),  it  follows  that 
this  coefficient  is  dimensionless.  Therefore, if  it  depends  upon  the  propagation 
constant  k,  it  must  be  a  f\mction  of  it  in  some  dimensionless  combination  ka  where 
a  has  the  dimension  of  length.  Let  us  choose  a  to  be  the  radius  of  curvature 
of  the  object  at  Q^  in  the  plane  of  incidence,  i.e.,  in  the  plane  containing  the 
incident  ray  and  the  normal  to  the  surface.  We  now  ass\ime  that  the  main  or  leading 
tern  in  D  (Q.)  does  not  depend  upon  any  other  geometrical  property  of  the  object. 
In  particular  it  does  not  depend  upon  the  radius  of  curvature  of  the  surface  in 
the  plane  normal  to  the  plane  of  incidence.  This  assumption  appears  to  be  in 
accord  with  the  ray  nature  of  the  field  propagation.  Of  course  small  additional 
terras  may  involve  other  geometrical  properties. 

The  decay  exponent  a  has  the  dimension  of  (length)   ,  We  will  assume  that  it 
depends  upon  k  and  upon  a,  the  radius  of  curvature  of  the  surface  in  the  normal 
plane  tangent  to  the  ray.  This  radius  is  the  reciprocal  of  the  geodesic  curvature 
of  the  ray.  We  will  see  in  a  number  of  examples  that  the  most  important  (i,e,, 
largest)  term  in  the  expression  for  a  involves  no  other  geometrical  property  of 
the  surface.  However , small  correction  terms  do  involve  other  geometrical  properties. 


-li- 
on the  basis  of  the  preceding  considerations  it  follows  that  the  main  terms 

in  D     and  a     can  be  detemdned  from  the  field  diffracted  by  any  object  of  simple 
m  m 

shape.     For  this  purpose  one  may  consider  diffraction  of  the  field  due  to  a  line 

source  by  a  circular  cylinder  of  radius  a.     If  the  line  source  is  parallel  to  the 

axis  of  the  cylinder  this  problem  simplifies  to  a  two-dimensional  one.     When  this 

field  is  compared  with  that  given  by  (12)   the   coefficients  D     and  a     can  be  found 

mm 

as  functions  of  a* 

In  order  to  effect  this  comparison  it  is  first  necessary  to  work  out  explicitly 
the  field  given  by  (12),     This  will  be  done  in  the  remainder  of  this  section.     Then 
it  is  necessary  to  obtain  another  expression  for  the  field  by  formulating  and 
solving  an  appropriate  boundary  value  problem  and  expanding  the  solution  asympto- 
tically for  large  values  of  ka.     This  is  done  in  Part  II,  Section  2,  for  fields 
satisfying  the  reduced  wave  equation  and  four  different  boundary  conditions. 
At  the  end  of  the  present  section  the  two  expressions  are   compared  and  a  table 

of  values  of  D     and  a     for  the  four  cases  is  obtained, 
m  m 

Consider  a  circular  cylinder  with  its  generators  parallel  to  the  z  axis. 
Let  an  (r,9)  coordinate   system  be  located  in  the  (x,y)  plane  with  0  *  0,  r  ■  p 
corresponding  to  the  location  of  a  line  source  (see  Figure  ii).     The  wavefronts 
associated  with  the  line  source  are  concentric  cylinders,  while  the  rays  are  the 
straierht  lines  emanating  from  the  source  and  orthogonal  to  the  wavefronts.     By 
using  these  rays  and  wavefronts  in  (2)  we  obtain  an  expression  for  the  incident 
field  produced  by  the  source.     If  R  is  the  perpendicular  distance  from  any  point 
to  the  source,   then  this  incident  field  at  the  point  is 

^^^'  "i    ■  H    ^I^^       * 

IXie  to  the  symmetry  of  the  problem,  it  is  clear  that  the  geodesies  are  arcs 

The  constant  in  (13)  has  been  adjusted  so  that  u.  represents  the  field  of  a  source 

■*■       i  (l) 
of  unit  strength,  i.e.,  one  for  which  the  exact  field  is  t-  H^   (kR), 
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of  the  generating  circle  of  the  cylinder.  Furthermore,  the  quantity  d<^Q-,)/(i6' 

has  the  value  unity  since  a  band  of  diffracted  rays  does  not  spread  on  the  cylinder. 

Also  the  radius  of  curvature  p,  in  (11)  is  infinite  so  that  p,/(p,+s)s  '  assumes 

-1/2 
its  limiting  value  s  '  .  Moreover,  because  the  geodesic  curvature  of  a  ray  is 

constant,  the  quantities  a  (t)  are  constants  and  D  (Q, )  -  D  (P,).  We  now  introduce 

m  m  1     m  1 

all  of  these  simplifications  into  (11)  and  use  (13)  for  u ,  (Q, )  noting  that  R  ?.t  Q,  is 
(p^-a^)^/^.  Then  (11)  be< 


)  comes 


-1/2 
(lU)    u^(r,e)  -  r8nks(p^-a2)^/2l    exp^k\^{p^'a^)^^^*i]^   ^1  ^  d; 


-1/2   r  .   o  o  wo  .  ..-1  _  2  ^^^-%^* 

'  e       , 
m 


Equation  (lU)  gives  the  field  associated  with  any  ray  from  Q  to  P  having  an 

arc  of  length  t  on  the  cylinder.  For  the  ray  QQ^P^P^the  length  t  has  the  value 

-1  -1  2  2  1/2 

t  ■  a©  -  a  cos  (a/p)  -  a  cos  (a/r),  while  s  »  (r  -a  )   ,  In  addition  to  this 

ray,  all  of  the  rays  which  are  tangent  to  the  cylinder  at  Q,  and  encircle  the  cylinder 

n  times  before  emerging  at  P,  are  also  diffracted  rays  through  P.  For  these  rays 

t  has  the  value  t  ■*-„■♦•  2nna.  Inserting  this  value  of  t  into  (lU)  and  summing 

over  n,we  find  that  the  diffracted  field  at  P  due  to  those  rays  which  are  tanpent 

at  Q,  is  given  by 

1/2    .  ,1/2, 


(15) 


f     r,\       /Q  ^x-1/2/  2     2^-1/U,   2     2x-l/U       Hu  /^r  2     2-,-^^^      r2     2t^'^\       in"! 
u^(r,e)  -  (8iik)     Mr  -a  )     ''*(p  -a  )     '^exp  ik  M P  "*  J       +    L^  -a  J       )  *  ft  \ 

5     (ik-a„)t^    r-  -1-1 

^     D^  e  "     °     1  -  exp[2n(ika-aaj^)]  . 

JB  L.  ^ 


Along  with  the  family  of  rays  described  above  we  also  have  rays  which  are 
tangent  to  the  cylinder  at  Qp,  encircle  the  cylinder  n  times  and  then  emerge  tan- 
gentially  at  P„  and  pass  through  P.  Computing  the  field  along  these  rays  yields 
a  result  of  the  form  (15)  with  ©  replaced  by  2r(  -  9,  Adding  this  result  to  (l5), 
we  obtain  for  the  total  diffracted  field  at  P 
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(     o^       /fi  u^-l/2/   2     2x-l/U,  2     2.-1/U        fT.    /r  2     2tV2      .2     2nV2N       i„  "l 
u^(r,e)   -  (8nk)     '    (r  -a  )     '^(p  -a  )     '^exp  ik  ([p  -a  J       +    [r  -a  J        J  ♦  ^ 


(16) 


'—      m 


expjjika-aa   )0j   ♦  exp[(ika-aa   )(2n-9)J 


m 


ipnCika-aajjj)! 


1  -  exp|2n(ika-aa   ) 


•  exp  -(ika-aajj|)(cos"  (^)  +  cos"  (^)) 


In  Part  II,  Section  2  wa  have  obtained  the  exact  solution  of  the  bovindary 

value  problem  for  diffraction  of  a  scalar  field  by  a  circular  cylinder.  The 

field  satisfies  the  reduced  wave  equation  and  satisfies  one  of  four  different 

conditions  at  the  cylinder  surface,  corresponding  to  different  physical  situations. 

We  have  also  expanded  these  solutions  asymptotically  for  ka  large.  In  all  four 

cases  the  leading  terns  in  the  expansion  of  the  diffracted  field  have  exactly 

the  form  piven  by  (16),  From  this  comparison  we  obtain  the  following  values  for 

a  and  D  in  these  cases 
n     m 


(17) 


m 


ik  -  iv  a 

m 


(18) 


D 


m 


■  e 


5ni/8  (2njlA 


■  AH^^^(ka) 
n 


-lV2 


The  expressions  on  the  right  sides  of  these  equations  are  defined  in  Part  IT, 

Section  2,  Values  of  a  and  D  ,  based  on  these  equations,  are  given  in  Table  I 

ID     m  f  IT. 

for  the  four  types  of  boundary  conditions  which  are  also  listed  in  the  table. 


5.  Electromagnetic  djf fraction 

In  the  electromagnetic  case  the  field  is  represented  by  two  vector  functions, 
the  electric  and  magnetic  fields.  Therefore,  as  we  pointed  out  above,  each  of 
these  fields  is  characterized  by  a  phase  ^^(s)  and  a  vector  amplitude  A(s).  We 
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assume  that  the  function  4  is  the  same  for  both  fields  and  is  given  by  the  method 

previously  described.  We  also  assvme  that  the  direction  of  each  amplitude  vector 

remains  constant  along  a  straight  ray  and  that  its  magnitude  is  given  by  the  previous 

construction.  Therefore, on  a  straight  ray  each  field  is  given  by  (2)  with  A  a 

vector.  Finally  we  assune  that  the  field  vectors  are  orthogonal  to  each  other 

and  to  the  ray. 

In  order  to  describe  the  field  on  a  diffracted  ray,  we  first  introduce  d,  i 

the  unit  tanpent  to  the  ray,  dp  the  outward  unit  normal  to  the  surface  and 

^-  "  t,   X  rtp.  Since  the  diffracted  ray  is  a  geodesic,  d^  lies  along  the  principal 

normal  to  the  ray  and  d.  lies  along  the  bn normal,  ftt  a  point  P  on  the  straifiht 

diffracted  ray  which  leaves  the  surface  at  P^  we  define  the  triad  of  vectors 

by  translating  those  at  P,  to  P.  We  denote  the  components  of  the  electric 

field  in  the  d,,  ^  and  d-  directions  by  E,,  E.  and  E,  and  use  a  similar 

notation  for  the  magnetic  field  H. 

We  now  assume  that  the  components  Ep  and  E^  propagate  independently  of  each 

other  and  obey  the  principles  satisfied  by  scalar  fields,  while  E,  is  zero.  There- 
fore the  diffracted  fields  E„  and  E-  are  each  given  by  (11)  with  A  exp  ik0.(Q-) 
replaced  by  E-CQ,)  or  E-(Q-).  However, the  diffraction  coeffidents  D  and  the 
decay  exponents  a  are  different  for  the  two  components. 

On  the  basis  of  the  above  assumptions  it  is  possible  to  express  the  diffracted 
electromagnetic  field  in  terns  of  two  scalar  fields.  First  we  determine  all  the 
diffracted  rays  which  pass  through  a  given  point  P,  Each  of  these  rays  originates 
at  a  point  0  on  the  surface  and  leaves  the  surface  at  some  point  P  ,  Then  we 
define  w(P,0  )  and  v(P,Q.)  as  the  values  at  P  of  two  scalar  fields  due  to  all 
diffracted  rays  which  orieinate  at  Q   and  leave  the  surface  at  P  ,  The 

corresponding  incident  fields  must  have  zero  phase  and  unit  amplitude  at  0  .  The 

field  V  is  constructed  by  using  the  values  of  D  and  a  appropriate  to  the  normal 

IB      m 

component  E^  while  w  is  constructed  with  the  values  of  D  and  a  appropriate  to 


the  taapential  component  E-, 
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Now  we  can  write  the  diffracted  field  as 

(19)       E^(P)  .  ^  [eJ(Q^)  v(P,0„)  d^CP^)  .  E^(Q^)  w(P,Q„)  ^3(P^)]  • 

In  (19)  Ej  and  E^  denote  the  components  of  the  incident  field  at  0  .  Equation 

(19)  is  just  an  analytic  statement  of  the  assunptions  described  above. 

For  a  perfectly  conducting  object  we  assume  that  D  and  a  are  the  same  for 

in     in 

the  tangential  component  E^  as  for  a  scalar  field  satisfying  the  condition  u  -  0 
on  the  surface.  These  values  are  given  in  Table  I,  Case  I.  Similarly, we  assume 
that  for  the  normal  component  Eg  they  are  the  saiae   as  for  a  scalar  field  satisfy- 
ing the  condition  dni/dn   -  0  on  the  surface.  These  values  are  also  given  in 
Table  I,  Case  II, 

For  an  imperfectly  conducting  object  we  assume  that  the  D^^  and  a^^   depend 
upon  the  torsion  as  well  as  the  radius  of  curvature  of  the  diffracted  ray.  If 
we  denote  the  torsion  of  the  diffracted  ray  at  the  surface  of  the  body  by  Z 

then  the  values  of  D  and  a  for  the  tangential  component  of  the  field  are  given 

m     m 

in  Table  I,  Case  IV  with 

bA^  '  \Ji  \i^     f(k,k^,at). 

Similarly,  the  values  of  a  and  D  for  the  normal  component  of  the  field  are  given 

'' '  mm 

in  Table  I,  Case  IV  with 

b/b^  -  l<^^ilk^^^^"■^f(k,k^,a'r). 

In  the  above  the  function  f  is  given  by  ^ 


f(k,k^,ar)     -     kN(k5  -k^)(ar)^  *  k|  I 


Notice  that  for    t  =  0,  i.e.,   for  plane  diffracted  rays,  the  function     f     is 
identically  equal  to  one.     Also, if  we  use  the  above   values  for  b^^  i"  Table  I 
and  then  let  k.   tend  to  infinity,we  see  that  we  obtain  our  previous  results  for 
perfectly  conducting  objects. 
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6,       Scalar  diffraction  by  a  sphere 

Let  UE  consider  the  diffraction  of  a  scalar  spherical  wave  by  a  sphere  of 
radius  a.     Let  the  sovirce  of  the  wave  beatr"p>a,  9«0ina  polar  coordinate 

system  with  the  origin  at  the  center  of  the  sphere   (see  Figure     5     )•     The  incident 

ikR 
field  u,    at  a  distance  R  from  the  source  is  found  from  (2)   to  be  e       /UnR.     The 

constant  has  been  chosen  to  correspond  to  a  source  of  unit  strength.     At  the 

2     2  1/2 
point  of  tangency  of  an  incident  ray  on  the  sphere,  R  ■   (p  -a  )  '    . 

To  apply  (ll)  we   observe  that,    just  as  in  the  case  of  the  cylinder,  the 

function  D     and  a     are  constants.     The  radius     a     entering  into  the  determination 
mm  "^ 

of  these  quantities  is  simply  the  radius     a     of  the  sphere.     The  ratio  d6"(Q- )/dcr(P- ) 
reduces  to  the  ratio  a  /a^   of  the  radii  of  the  latitude  circles  at  points  Q^   and  P, 
respectively.     In  order  to  facilitate  the  calculation  of  this  quantity  we  have 
drawn  in  Figure   5  a  cross  section  of  the  sphere  through  the  plane  containing 
the  diffracted  rayr.    A  simple  calculation  then  shows  that 

(20)  a^  -  a(p2-a2)^/2  p-1 

(21)  a,   -  ar"  (a  sin  9  -  s  cos  9)  . 

From  Figure  5  it  is  also  apparent  that  the  distance  t  which  is  the  length  of  the 
arc  of  the  circle  from  Q,  to  P.  is  given  by 

(22)  t  ■  a©  -  a  cos"  (a/r)  -  a  cos~  (a/p)  , 

In  Appendix  I   it  is  shown  that  p.,  the  non-zero  principal  radius  of  curvature  of 
the  diffracted  wavefront  on  the  body,  is 

(23)  p,     ■     a(a  sin  9  -  s  cos  9)(a  cos  9  +  s  sin  9)      . 
From  (20),    (21),   and   (23)  we  obtain 


Figure?:     The  two  families 
of  rays  emanating  from  a  point 
source  located   at  Q  and  diffracted 
around  a  sphere  to  a  point  P  in 
the  dark  region.     The   cross-section 
of  the  sphere   shown  is  the  plane  of 
the   polar  axis  and   the  line   Joining 
Q  and  P.     One   family  of  rays  is  in- 
cident at  Q-    and  leaves  at  P, ,  while 
the  other   family  is  incident  at  Q« 
and   leaves  at  P„ . 


^^g^re  6t     The   two  families  of  rays 
emanating  fron  a  pcint  source  located 
at  Q  and  diffracted  around  a  sphere  to 
a  point  P  in  the  illuminated  region. 
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Upon  substituting  (22)  and  (2h)  together  with  the  value  of  the  incident  field  on 
the  body  into  (11)  we  obtain 


exp  (ik-a  )a  |9  -  cos  (a/r)-  cos"  (a/p)| 


To  obtain  the  complete  expression  for  the  diffracted  field  at  P,  we  must 
include  diffracted  rays  which  are  tangent  at  Q^  and  encircle  the  sphere  n  times 
before  leaving  at  P,,  Then  we  must  also  include  a  similar  family  of  rays  tangent 
at  Q^  leaving  the  sphere  at  Pp.  As  in  the  case  of  the  cylinder,  the  field  at  P 
due  to  all  the  rays  tangent  at  Q^  can  be  obtained  by  summing  a  geometrical  series, 
In  the  case  of  a  sphere,  however,  a  phenomenon  which  is  not  present  in  the  case 
of  the  cylinder  appears.  All  of  the  diffracted  rays,  except  the  ones  sketched, 
pass  through  foci  at  the  north  and  south  poles.  We  postulate  that  the  phase  de- 
creases by  ii/2  as  a  ray  passes  through  a  focus,  in  accordance  with  optical  experi- 
ence. The  series  which  must  be  summed  is  then  exp  2nna( ik-a  )-inn  from  n  -  0 
to  n  «  00.  Upon  summing  this  series  and  using  the  result  in  (25),  we  obtain  for 
the  contribution  to  the  field  at  P  from  all  these  rays  which  are  tangent  at  Q, 
and  encircle  the  sphere  any  number  of  times 

''"     "  ■  U(pr|.ioe|)V^(r^-a^)V'.(p^.a^)lA  f  ^j^  *  exp(2,a(ik-a„)J 


exp  (ik-a  )a  ■j  e  -  cos'  (a/r)  -  cos  (a/p)i 


The  field  due  to  rays  *«ngent  at  Q2  is  obtained  from  (26)  by  replacing  9  by  2n-e 
and  multiplying  the  resulting  expression  by  -i.  Thus  the  total  diffracted  field 
is 
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(27)  u^(r,e)  -  w(r,9)  -  i  w(r,2n-©). 

I, 

When  use  is  made  of  (1?)  and  (13),  it  is  easily  seen  that  (27)  agrees  exactly 
with  (10)  of  Part  II  which  was  obtained  by  asymptotically  expanding  the 
solution  of  the  corresponding  boundary  value  problem. 

We  have  tacitly  assumed  (see  Figure  5),  in  constructing  the  diffracted  field 
as  given  by  (27), that  the  point  of  observation  is  in  the  geometrically  dark  region. 
As  a  further  check  on  the  geometrical  construction  of  the  diffracted  field  we  con- 
struct the  diffracted  field  in  the  geometrically  lit  region.  The  geometry  is 
shown  in  Figure  6, 

Simple  geometric  calculations  show  that  in  this  case  the  distance  t  from 
Q.  to  P-  along  the  circle  has  the  value 

t   ■  (2n-9)  a  -  a  cos  (a/r)  -  a  cos"  (a/p). 

The  value  of  a  is  again  given  by  (20),  while 

(28)  a,   -  ar*  (a  sin  0  +  s  cos  9), 

In  Appendix  I  it  is  shown  that  p,  is  equal  to  the  stradght  line  distance  between 
P-  and  C, .  From  Figure  6  this  is  given  by 

(29)  p.  -  a(a  sin  9  +  s  cos  e)(s  sin  9  -  a  cos  9)   • 

When  we  substitute  (20),  (28)  and  (29)  into  the  expression  for  the  amplitude 
it  becomes  identical  with  (2h). 

As  before, rays  which  are  incident  at  Q^  and  encircle  the  sphere  n  times  before 
emerging  at  P  are  also  diffracted  rays.  Since  each  of  these  rays  passes  through  an 

odd  number  of  foci  the  phase  at  P  along  the  n   number  of  this  family  is 


(30) 


iki(p^-a^)^^^  *  (r^-a^)^/4*  (ik-Q^)|t^  *  2(n-l)7ial  -  i(2n+l)n/2 


Upon  using  (30),    (2h)   and  A     ■   (Un)"     in   (11)   and   summing  over  n  from  n  =  1  to   oo 
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we  obtain  as  the  contribution  to  the  field  at  P  from  all  these  rays  which  are 
tanpent  at  Q, 

(31)  v(r,9) -,      ,    .^,,1/;,   ;     ;,l/l.,   ;     i,\/k      E  »J'-  ♦  «xp|2«a(lk-a  )J 

Un(pr|sin6|)   '    (r  -a  )   '    (p  -a  )   '  m  l-  -J 


exp  (ik-a   )a  ^2n-e-cos'   (a/r)-cos"'^(a/p)  > 


Since  the  diffracted  rays  through  P  vrtiich  are  initiated  at  Qp  pass  through 
a  focus  off  the  body  (C, ),we  must  also  include  a  phase  decrease  of  v/2  when  passing 
through  this  focus.     Thus  each  member  of  this  family  of  rays  undergoes  a  phase  de- 
crease which  is  an  even  multiple  of  n.     Upon  takinp  this  phase  shift  into  consider- 
ation we  find  that  the  field  due  to  the  rays  tangent  at  Q^  is  obtained  from  (31)  by 
replacing  0  by  -0  and  multiplying  the  resulting  expression  by  -i.     Thus   the  total 
diffracted   field  is 
(32)  u^(r,0)     -     v(r,e)  -  i  v(r,-0). 

Again^by  using  (1?)  and  (18),  we  see  that  (32)  agrees  with  the  diffracted  field 
in  (11)  of  Part  II  which  was  obtained  by  asymptotically  expanding  the  solution 
of  the  corresponding  boundary  value  problem. 

7.       Scalar  diffraction  by  an  arbitrary  cylinder 

In  this  section  we  treat  scailar  diffraction  of  an  obliquely  incident  plane 
wave  by  a  cylinder  of  arbitrary  convex  cross  section.     We  will  find  that  if  the 
geometrical  theory  is  valid  for  a  normally  incident  plane  wave,   it  is  also  valid 
for  an  obliquely  incident  one.     To  do  this  we  first  indicate   the  manner  in  which 
the   corresponding  boundary  value  problem  at  oblique  incidence  reduces  to  one  at 
normal  incidence.     Then  we  show  that  the  solution  obtained  by  the  geometrical 
theory  reduces  in  exactly  the  same  way. 

First  we  introduce  a  rectangular  coordinate   system  with  the   z-axis  parallel 
to  the  generators  of  the  cylinder .     We  orient  the  x-axis  so  that  the  incident 
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plane  wave,  assumed  to  be  of  unit  amplitude,  is  given  by 

„     «ik(z  cos  a  *  J   sin  o) 
U.  "  e  "'       • 

We  also  introduce  a  curvilirtear  orthogonal  coordinate  system  (h,,h^)  in  the   (x,y) 
plane  with  h-    ■  a     being  the  equation  of  the  generating  cuirve  of  the  cylinder. 
Let  us  suppose  that  U  is  a  solution  of  the  following  equations 

(V^  +  k^)U     -     0,  h^     >     a 

(\7^  +  k^)U     -     0,  h^    <     a. 

Furthermore  we  suppose  that  at  h.    -  a,  U   satisfies  one  of  the  boundary  conditions 
listed  in  Table  I,     The  incident  field  U.   is  given  above  and  the  scattered  field 
U-U.   must  satisfy  the  appropriate  two-dimensional  radiation  condition.     Now  we 
define  V(h^,h2)  by  setting  U(h^,h2,z)  -  V(h^,h2)e^^^  ^®°.     Then  it  follows  that 
V  satisfies 


(33) 


(V^  +  k'^)V     =     0,  h^     >     a 

(7^  +  k^^)V     -     0,  h^     <     a. 

It  also  follows  that  V  satisfies  the  same  boundary  conditions  as  does  U  at  h^  ■  a, 

ik'v 
that  V,  ■  e   "^  and  that  V-V.  satisfies  the  two-dimensional  radiation  condition. 

t     2    2    2  1/2 
In  (33)  k'  -  k  sin  a  and  k,  -  (k,  -  k  cos  a)  '  ,  while  7,  is  the  two-dimensional 

Laplacian  in  the  variables  h,  and  h^.  These  equations  show  that  V  is  the  two- 

ik'y 
dimensional  field  in  a  medium  of  propagation  constant  k'  due  to  the  wave  e 


incident  on  a  cylinder  having  propagation  constant  k- .  From  it  we  can  find  the 
diffracted  field  U  by  multiplying  V  by  e        and  replacing  k'  by  k  sin  a  and 
kJ  by  (kf  -  k  cos  a)  '  .  We  will  now  construct  u  and  v,  the  fields  given  in  these 


Care  must  be  exercised,  however,  in  Case  III.  Because   k  appears  in  the 
boundary  condition,  we  must  also  replace  Z  by  Z/sin  a  in  Table  I, 
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two  problems  by  the  geometrical  theory.  We  will  then  show  that  u  can  be  obtained 
from  V  in  the  same  manner  as  U  can  be  obtained  from  V,  This  demonstrates  the 

consistency  of  the  geometrical  theory.  It  also  permits  the  solution,  by  the  peo- 
metrical  theory,  of  a  three-dimensional  problem  to  be  reduced  to  that  of  a  two- 
dimensional  one. 

To  construct  u  let  the  cylinder  be  given  by 

(3U)         x(t,z)  -  (x(t),  y(t),  z). 

Here  x(t)  and  y(t)  are  such  that  h- ^x(t),y(t))  -  a,  and  t  is  arc  lenp^th  alonp  this 
curve.  Then  the  geodesies  on  the  cylinder  are  helices  given  by 

(35)  x(s)  -  (xfps+3),  y(ps+p),  YS+6). 
The  arc  length  along  each  helix  is  s  if 

(36)  P^  +  Y^  -  1. 

The  constants  p,  p,  y,  and  5  depend  upon  the  position  and  direction  of  the  geodesic 
at  E  "  0,  To  evaluate  them  we  let  (x,y,z)  be  the  coordinates  of  the  point  of  ob- 
servation and  impose  the  condition  that  at  s  ■  0  the  tangent  to  the  helix  points 
in  the  direction  opposite  to  that  through  (x,y,z).  This  condition  yields 

(37)  xO)  -  h  p  M?)     "     X 

(38)  y(p)  -  h  p  y(p)  -  y 

(39)  5  -  ..  Y  .  z. 

In  these  equations  h  is  a  constant  of  proportionality  determining  the  distance 
along  the  tangent  to  the  point  of  observation  and  the  dot  represents  differentiation 

with  respect  to  the  argument. 

•J 

We  are  assuming  that  the  helix  is  parametrized  so  that  the  point  s  ■  0  is  the  one 

at  which  the  diffracted  ray  leaves  the  cylinder  and  proceeds  to  the  point  of  ob- 
servation. As  s  increases, the  corresponding  point  moves  backwards  along  the 
diffracted  ray  to  its  point  of  initiation. 
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We  now  denote  by  s     the  fyeodesic  distance  along  the  cylinder  from  the  point 
of  incidence  to  the  point  at  which  the  diffracted  ray  leaves  the  surface.     Then 
at  s  »  s     the  tangent  to  the  helix  must  point  in  the  direction  opposite  to  the 
direction  of  incidence.     This  requires  that 

(iiO)  p  x(p  s^  +  p)     -     0 

(ill)  p  y(p  s^  +  P)     =•     -  sin  a 


(U2)  Y     ■     "  cos  a. 

Since  the  argument  of  x  and  y  is  the  arc  length  of  the  generating  curve  of 
th<?  cylinder  we  have  x  +  y  -  1,  Therefore, upon  squaring  and  adding  (37) 
and  (38)  we  obtain 
(U3)         h  «  rpipf-^  . 

Here  we  have  set  r.  »  (x(p)-xy  +  (y(p)-yy   j   .  To  analyze  (Uo)  we  observe 
that  the  convexity  of  the  generating  curve  ensures  that  there  are  precisely  two 
points  on  the  curve  for  which  x(t)  "0.  Therefore,  if  t,  and  tp  are  two  distinct 
points  at  which  x(t)  ■  0,  and  T  denotes  the  length  of  the  generating  curve,  (IjD) 
permits  us  to  write 


•i" 


ihh)  p  sj;-"^  +  p  -  t^  +  nT 


(U$)         P  s^^^  +  p  -  tg  +  nT. 


•2       .2 


From  the  fact  that  x'  +  y'  -  1  we  have  y  ■  +  1  at  t^  and  tj.  We  let  t^  be  such 

that  y(t-)  -  1  and  then  y(t2)  -  -1.  Then  from  (hi),  the  following  two  values  of 

p  correspond  to  s    and  s^   respectively 
n       n 

(U6)         p  *     -   sin  a 

(2) 
(li7)  p     «  sin  a  » 
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Upon  substituting  (U6)  or  (U7)  into  (U3)  we  obtain 

(U8)  h     -     r-  /  sin  a  . 

P 

When  this  value  of  h  and  the  above  values  of  p  are  inserted  into  (37)  and  (38), 
these  equations  become 

(U9)         x(p)  ±  Tp  x(?)     -  X 
(5t))         y(p)  ±  Tp  y(p)  =  y. 

These  equations  have  two  roots  p,  and  pp*  corresponding  to  the  upper  and  lower 
signs  respectively  .  The  value  of  y  is  given  by  (U2)  while  from  (39) 

(51)  6  ■  z  -  Tq  cot  a  , 

P 

Upon  substituting  the  value  of  p^       into  (UU)  and  the  value  of  p^       into  (,hS)  vre 
obtain  for  each  value  of  n, 

(52)  s^-*-^     -     (p^  ■  H  *  "'^^   /  ^^"  ° 

(53)  s^^^     -     (-pg  *  tj  +  nT)  /  sin  a. 

The  distance  s  appearing  in  (11)  has  two  values,  one  for  each  family  of  rays, 
i.e.,  for  each  value  of  p,  A  simple  calculation  shows  that  these  values  are  given 
by 
(5U)         s^  ■  r^  /  sin  a  i  «  1,  2, 

It  is  also  easily  seen  that  a  (s),  the  radius  of  curvature  of  the  helix,  is  given 
in  terms  of  a  (t),  the  radius  of  curvature  of  the  generating  carve,  by 

As  a  consequence  the  value  of  rp  depends  upon  the  value  of  p.  For  p  -  p,  we 

denote  r„  by  r^ ,  while  for  p  ■  p_  we  denote  r^  by  r„ .  When  no  confusion  can 
p     J.  ^  p  ^ 

occur  we  will  continue  to  use  r_» 

p 
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(5$)  a  (s)  -  -Jy-  a  (ps  +  p). 

^       sin  a   ^ 

To  complete  the  calculation  of  the  geometrical  quantities  we  need  the  value 

of  the  incident  field  at  the  points  Q  where  the  diffracted  rays  are  initiated 

a- 
on  the  surface  •  The  values  of  U.  at  the  points  Q  are  determined  by  substitutino; 

the  values  of  y  and  z  at  the  points  of  incidence  into  e^^^^"^^"®  *  zcosa)^  Making 

use  of  our  previous  calculation  we  obtain 

(56)     y(i^')     -  y(p^^^^n^^  *Pi)  -  y(ti-nT).  y(t.) 


(57) 


(•^n^O  '  ^^n^^  *  ^^^^  "  ^"  [""i  *  "^  *  ^"^^^^i  '  (-l)^Pi]'=°^  ° 


Because  of  the  cylindrical  sjnrmetry  the  quantity  dotQ)/datP)  in  (11)  is 
equal  to  1,  while  p,  ■oo.  Consequently,  taking  into  account  all  of  the  rays 
through  P,   (11)  becomes 

m  n«C  •-  J  ~J 


«.(^2)^i'^'  E„  ^„(«^'>i€')=^f  ^(4''".)-  j  "%(='^ 


o 

'    '  "si 


'o 

Before  substituting  the  values  calculated  for  the  quantities  appearing  in  (58) 
we  observe  that  D  (O^'''  )  is  independent  of  n.  This  is  because  it  depends  only 
upon  a  (s   )  which  is  clearly  independent  of  n  as  we   see   front  (55),   (UU)   and   (U5). 

■? ~~~  ~ 

There  are  two  families  of  diffracted  rays,  cne  correspondinp  to  each  value  of  p. 

The  rays  of  one  family  are  incident  at  the  points  Q^   defined  by  (lih).  They  all 

leave  the  cylinder  at  the  point  P,  defined  by  (U9)  and  (53)  with  P  *  p, .  The 

(2) 
rays  of  the  other  f airily  are  incident  at  the  points  Q^   defined  by  (U5).  They 

all  leave  the  cylinder  at  the  point  Pg  defined  by  (U9)  and  (50)  with  P  -  Pg. 
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Similarly,  the  function  a  (s)  which  also  implicitly  deperds  upon  the  points 

g^^^  at  which  the  rays  are  initiated  is  also  independent  of  n.  Carrying  out 
n 

the  substitutions  and  summing  over  n  we  obtain 

1 

(59)      u^(F)  -  (sina)^  e^^^^^'^  JZ  \^hK(S^^^^>l^^ 

n 

(p^-t^)/sina 

exp  iksina  (p^-t^+r^^+yCt^))  -     aj^(s,Q'-'-^^ds 

i 
T/sina 


1  -  expfikTsina  -    a^(s,Q^-'-^dsl 


o  ^ 

In  this  result  the  field  associated  with  the  rays  whidi  leave  the  surface 
at  P,   has  been  written  explicitly.     The    *  denotes  the  corresponding  field  on  the 
rays  which  leave  the  surface  at  P-,     The  expression  for  this  field  can  be  obtained 
from  the  above  sum  by  replacing  the  subscripts  and  superscripts  1  by  2  and 
chanping  the  signs  of  p.    and  t^  except  in  the  argument  of  yCt^),     The  syrnbol 
a  (s,Q      )   denotes  the  value  of  a     at  a  point  a  distance  s  along  a  ray  from 
any  of  the  points  Q,         at  which  the  rays  are  initiated.     Similarly,   the  symbol 

D  (Q      )  denotes  the  value  of  D     at  any  of  the  points  Q/      .     The  values  of  a 
m^-'  mn  n 

and  D     can  be  obtained  from  Table  I, 
m 

Equation  (59)  is  the  diffracted  field  due  to  oblique  incidence  of  a  plane 
wave  on  a  cylinder  of  arbitrary  crofs  section.     If  we  set  a  -  n/2  in  this  equation, 
WB  obtain  the  result  for  normal  incidence  previously  obtained  in    [2] .     Now  if  we 

replace  k  by  k     and  k_   by  k^  in  the  result  for  normal  incidence,  we  obtain  the 

t 

function  v.  We  now  wish  to  show  that  from  v  we  can  obtain  u  by  replacing  k  by 

k  sin  a,  k^  by  (k-  -  k  cos  a)  ^     and  multiplying  the  result  by  e^  ^'^  °,  Upon 
performing  these  operations  in  (59), we  find  that  u  and  v  are  related  as  above 
provided  that 


27  - 


W/sina  ^         /  1/2 

(60)  a  (sjk,lc,)ds     »  a    (  sjk  sina,    \^-,~^  cos  a  ]  ds 


(61)  (sin  a)-^/*^  Dj^(A;k,k:^)     -     T>^  [ajK  sina,    fk^-k^oos^a]       )      . 

In  (61)  A  denotes  an  arbitrary  point  while  in  (60)  W  denotes  an  arbitrary  upper 

llndt  of  integration.  The  dependence  of  D  and  of  a  on  k  and  k,  has  been  in- 

m        m         J. 

dicated  explicitly. 

To  verify  (60)  we  express  the  value  of  a  ■  a  (s)  appearing  in  the  expression 
for  a  in  Table  I  in  terms  of  a  by  means  of  equation  (55)  to  obtain 

(62)  a  (sjk,k,)  .  e-i"/6  f -^ ^""^  a     si^^/^a  • 

Upon  substituting  (62)  into  the  left-hand  side  of  (60)   and  changing  variables 
by  setting  s  sin  a  "  E,  ve   obtain 

(63)  a Jsik,k,)ds     3     e-i"/6  f}^y  ^ . 

In  (63)  we  have  used  the  fact  that  p  ■  +  sin  a. 

Now  in  Cases  I   and  II,  where  a     is  a  numerical  constant,  equation   (63) 
demonstrates  the  validity  of  (60).     In  Cases  III  and   IV,  however,  q^  is  a  function 
of  k,  k-,  and  a.     Therefore,  the  right-hand  sides  of  (63)   and  (60)  will  be  identical 
if 

a/2 


(6b) 


/  r2     2       2l^^^         \ 

q^(k,k^,8   )   -  q^  I  k  sine,     k^-k  cos  a         >  *d  ) 


In  Case  Til,  when  q^  is  also  a  function  of  Z,  the  Z  contained  in  the  right  hand 
side  of  (6li)  must  be  replaced  by  Z/sina  (see  footnote, p,  21),  When  use  is  made 
of  (55),  together  with  the  entries  for  A  {q^)/k{q^)   in  Table  I,  (6U)  is  easily 
verified. 
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To  verify  (61)  we  first  use  Table  I  to  write  Djjj(A;k,k^)  in  the  form 


(65)  DjA}k,k^)     -     k"^/^a^/^  f(q^). 


In  (65),  ^(%.)   is  in  all  four  cases  a  function  only  of  the  variable  a  ■  a  (k,k,,a  ), 
When  use  is  made  of  (55)  and  (6h),  (65)  is  readily  verified. 

Thus  we  have  shown  that  the  result  (59),  given  by  the  geometrical  theor>' 
for  diffraction  by  an  obliquely  incident  wave,  can  be  obtained  from  tho  result 
for  a  normally  incident  wave.  The  relation  between  these  two  fields  is  the 
same  as  for  the  solutions  of  the  corresponding  boundary  vaD.ue  problems  con- 
sidered at  the  beginning  of  this  section, 
8,   Electromagnetic  diffraction  by  a  sphere 

Let  us  now  consider  the  diffraction  of  a  linearly  polarized  plane  electro- 
ma''netic  wave  by  a  sphere.  We  will  use  the  same  notation  as  in  Section  6  and 
assume  that  the  incident  electric  field  is  parallel  to  the  x-axis  and  propagating 
in  the  positive  z  direction.  Then  the  incident  electric  field  can  be  written 


as 


(66)  E^  -  Ae^^^  E^  .  E^  =  0, 

The  incident  rays  are  parallel  to  the  z-axis  so  the  only  points  of  tangency  will 

lie  on  the    "equator"    of  the  sphere  in  the  plane  z  »  0.     The  diffracted  rays 

are  plane  curves  (great  circles)  on  the  sphere.     Therefore, the  rays  through  a 

point  P  with  coordinates  (r,0,0)  will  lie  in  the  plaine  0  »  constant  whi.ch  passes 
through  P,     This  plane  is  shown  in  Figure  7,     There  are  two  families  of  diffracted 

rays  thro\;igh  P.     One  is  initiated  at  Q^(a,Ti/2,0)   and  leaves  the  sphere  at 

P-Qa,9+cos     (a/r),0)  provided  e+cos     (a/r)  <  n.     The  other  family  originates 

at  Q-(a,n/2,0+n)   and  leaves  the  sphere  at  Pp^a,9-cos~  (a/r),(i^  prov''^"'.- 

&  >  cos     (a/r).     If  9+cos~  (a/r)  >  n  then  P.    has  coordinates    (a,2n-^-cos~^(a/r)  jfi+n), 

If  9  <  cos"  (a/r)  then  P       has  coordinates  ^a,cos     (a/r)-e,^+r).     The  rays  of 

each   family  encircle  the   sphere  any  number  of  times  before  leaving  it. 

We  will  calculate  the  diffracted  field  E^(P)  by  using  (19).     Therefore, we 
first  calculate  the  rectangular  components  of  the  d  vectors  at  the  points  Q^  and  Q2' 


28«  - 


(rA^) 


Figure  7 1      Th#  ray  system  associated  with  the  diffraction 
of  an  electromagnetic  wat/e  by  a   sphere.     The 
diffracted  rays  are  incident  on  the  equator 
at  the  points  Q,    and  Q^  and   leave  tangontially 
at  the  points  P.   and  P^, 
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These  vectors  are 


(67  )   d]^(Q;^)  -  (0,0,1),  d2(Q^)  -  (co80,6in0,O),  "^AQj)   '   (-sin0,cos0,O) 

(68)    ^1(02^  "  (0,0,1),  d2(Q2)  =•  (-cos0,-sin0,O),  d^(Q2)  "  (sin0,-cos0,O) . 
Similarly,  the  spherical  components  of  the  corresponding  vectors  at  the  points  P, 


and  Pp  are 


d^(P^)  .  (r-l(r2-a2)^/^-a/r,0),  d^{?^)   -  (a/r ^HrW )^^\6) , 

^  (P  )  -  (0,0,1) 
(69) 

d^(P2)  -  (r-^(r2-a2)^/^a/r,0),  d2(P2)  -  (a/r,-r-^(r^-a^)^/^o), 

d3(P2)  -  (0,0,-1)  . 

From  (67)  and  (68)  and  the  expression  (66)  for  the  incident  field,  we  find  that 
at  Q^  and  Q^  this  field  has  the  components 

^2^\^     '     A  cos  0,      E2(Q2)  -  -A  cos  ft), 
(70) 

E^(Q^)  -  -A  sin  0,     ^\(^)     '     A  sin  0. 

By  using  the  methods  of  Section  6  it  can  easily  be  seen  that  the  function  v(P,Ot ) 
appearing  in  (19)  is  identical  with  the  function  v(r,9)  defined  by  (UO)  of  Part  II. 
Similarly,  v(P,Q2)  is  equal  to  e'^"'  v(r,-e);  while  w(P,Q^)  is  equal  to  the  function 
u(r,e)  also  defined  by  (liO)  and  w(P,Q2)  is  equal  to  e'^"'^u(r,-©). 
Upon  using  the  above  together  with  (70)  in  (19)  we  obtsdn 

(71)     E^(P)  -  A  cos  (jf[t^{?^)v{T,9)  '   6"^"/^  t^{?^Mr,'9)\ 

'  A  sin  0r?3(P^)u(r,e)  -   e"^"/^  d3(P2)u(r,-e)l    . 
If  we  now  express  the  d  vectors  in  terms  of  their  spherical  conponents  as  given 
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by  (69)  W8  see  that  (71)  is  in  agreement  with  (Ul),  ()i2),  and  (U5)  of  Part  II 
which  was  obtained  by  asymptotically  expanding  the  solution  to  the  corresponding 
boundary  value  problem. 

9.  Electromagnetic  diffraction  by  an  arbitrary  cylinder 

We  now  consider  the  diffraction  of  a  plane  linearly  polarized  electro- 
magnetic wave  by  a  cylinder  of  arbitrary  convex  cross  section.  We  introduce  a 
rectangular  coordinate  system  as  in  Section  7  so  that  the  direction  of  incidence 
is  given  by  the  irnit  vector  d^  ■  (0,  sin  a,  cos  a).  If  the  electric  vector 
makes  an  angle  (n/2  -  6)  with  the  (x,y)  plane, then  the  rectangular  components 
of  the  incident  electric  field  can  be  written  as 

E^  -  A  sin  e  eik(y8ina  ♦  zcosa) 

X 

/_„x  _i    ,  a       ik(ysina  +  zcosa) 

(72)  F  "  A  cos  B   cos  a  e   '' 

y 

E^  -  -A  cos  6  sin  a  e^^^^^*^  *  ^°°«°^  • 

t 

The  diffracted  ray  system  has  been  described  in  Section  7.  At  a  piven 
point  P  it  consists  of  two  families  of  rays.  The  rays  of  one  family  are 

initiated  at  points  q}       and  leave  the  cylinder  at  the  point  P.,  while  the  rays 

(2) 
of  the  other  family  are  initiated  at  points  Q    and  leave  the  cylinder  at  the 

point  Pp. 

In  order  to  calculate  the  diffracted  field  by  means  of  (19)  we  first 

calculate  the  rectangular  components  of  the  d  vectors  at  the  points  Q    and 

(2) 

Ql      *     These  vectors  are  independent  of  n  and  are  given  by 

i^(Qjj)  -  (0,  sina,  cosa),  ^^(qJ)  -  (1,0,0),   ^3(0^;)  -  (0,  cosa,  -sina) 

(73) 

^(Q^)  -  (0,  sina,  cose),  ^2^<^^  '   (-1>0>°)»  ^^^^n^  "   (°»  -cosa,  sina). 

Similarly,  in  the  notation  of  Section  7,  the  rectangular  components  of  the 
corresponding  vectors  at  the  points  P,  and  Pp  are 
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^(P^)  -   {^{^^)  Sim,  f(^^)  Sim,   coso),     ^(P^)  -  -&^ip^)^{?^),  y(p^),  o) 

3^(P^)  -   (ap(p^)f(p^)  cosa,   -ap(p^)  x  (p^^)  cosa,  -sina) 

(7U) 

ci^CPg)  -   ^xCPg)  sina,   -^(pg)   sina,    cosa),   t^{?^)  =  -ap(p2)(K  (P2^»  ^^^2^'   °) 


ctj(P2)  ■(ap(p2^  3'^P2^  °°^*  "^p^^2^  ^'^^2^  °°^°'    ^^"^V 


(1) 


By  using  (72)  and  (73)  we  find  that  the  components  of  E     at  the  points  Q 

(2) 
and  Q^   '^  are 
n 

EW'^h  -  0         ,  eJ(q(^^)  -  A  3in9  e^^^^^^-^^  *  ^  ^°^^^ 
In  *     2     n 

(75)  E^(Q^^^  -  A  cos9  e^^^^y^^"^  *  ^  """^^ 

3     n     ' 

E^Q^^)),  (3  gi(Q(2))  .  .^  g.^e  ^ik(y3ina  *  z  cosa) 

i     n  t     n 

(76)  Ei(Q^2)j  .  .^  ^^gQ  3ik(ysina  +  z  cosa)  . 

3     n 

In  (75),  y  and  z  are  to  be  evaluated  at  Q^     ,  while  in  (  76)  y  and  z  are  to 

(2) 
be  evaluated  at  Q^   \     The  values  of  yand  z  at  the  points  Q     are  given  by  ('56) 

and  (57)  of  Section  7. 

In  the  present  notation  (19)  becomes 

(77)  E^(P)  -y    ^     r4(^?^)  '  (P,Qp^)  t(Pj  +E^(Qp^)    w(P,Qp^(P  ~    » 

j^l  n^O  L*^     "^  n  <;     J  3     n  n        j     j  j 


The  calculations  of  Section  7  show  that 
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(78) 


n-0  '^ 


r.    4(Qi^^)w(P,Q^-^^)  -w^(P)  -  (-D^^^Acose 


y  (sina)  e 


?  _ikzcoso 


^D^(Pj)D^(Q^J^)r'^/'^exp  iksina  (rj+(-l)Jt^-(-l)^Pj^y(tj))- 


1 


fT-l)^t  -(-l)Jp/]/sina 


a^(s,Q^J^)d8 


(J). 


1  -  exp 


T/sina 


.(j)' 


ikTsino  -  |   Oj^CsjQ^^Ods 


-1 


The  values  of  D  and  a  in  the  expression  for  v  (F)  are  those  appropriate  for 
the  normal  component  of  the  electric  field,  while  the  values  of  these  quantities 
appearing  in  the  expression  for  w.(P)  are  those  appropriate  for  the  tangential 
component  of  the  electric  field.  In  Section  5  we  have  described  the  way  in  which 
quantities  are  to  be  determined. 

Upon  using  (7U),  (77),  and  (78)  wb  obtain  for  the  rectangular  compcnents 
of  the  diffracted  field 

E^(P)  -  JZ   [-v^(P)ap(Pj)  XO.)  *  w^(P)apOj)  ^(Pj)cos  aj 
(79)     fJ(F)  .  ^   r-v^(F)apOj)  y(|3j)  -  v.^(F)apOj)  iO^jcos  J 


<<^) 


■■^ 


(-1)^  Wj(P)sin  a   . 


Let  us  now  specialize  (79)  to  the  case  of  a  point  P(p,/&,z)  in  the 
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shadow  of  an  imperfectly  conducting  circular  cylinder  of  radius  a.     In  this 

case  the  functicn  x(t)  is  equal  to  a  cos(t/a),  while  y(t)  is  equal  to  a  sin(t/a). 

Simple  calculations  show  that  p,    ■  a(^-f)   and  ^^  '  a(0* '^)  where    f  ■  cos     (a/p), 

2     2  1/2 
and  0  <  t<  r/2.     The  quantity  r.  is  equal  to  (p  -a  )  '    ,  while  t,    ■  0  and  tp  ■  an, 

The  radius  of  curvature  a     is  a,  while  the  torsion     f   which  enters  into  the 

P 

determination  of  c  and  D  is  equal  to  a  sina  cosa. 
m      m 

In  order  to  identify  (79)  with  the  asymptotic  expansion  of  the  exact 

solution  we  first  determine  the  values  of  a  and  D  appearing  in  (78).  From  the 

m  n 

results  of  Section  5  we  see  that  for  the  function  v.(P)  the  values  <af  a     are 
determined  by  the  equation 


-      ^-.-(^^' 


A,x? 


Q)- 


T   kV 


1 


(kj-k^)(a^r)^  *  k^V      , 


where  we  have  denoted  the  radius  of  curvature  of  the  diffracted  ray  on  the 

surface  bv  a     in  order  to  avoid  confusion  with  the  radius  a  of  the  cylinder. 
? 

Upon  substituting  a(sin  a)       for  a     and  a     sina  cosa  for    T,  we  may  write   (80)   as 

(81)     "''"'  .  ,^(''"•^^^(''^f  ii^"''"! 


1 


2 


where  we  have  introduced  k'  ■  k  sin  a  and  k,  »  k^  -  k  cos  a.  From  (81)  we 


d. 


see  that  q     entering  into  the  determination  of  v, (P)  are  identical  to 
the  q^  entering  into  the  determination  of  the  function  v(p,0)  of  Part  II,  .Section 
U.  In  fact,  upon  ''oraparing  (78)  and  (2U)  of  Part  II  we  have 


(82) 


vJ(P)  -  A  sin  e  e^''^°°^°v(p,0) 


d/r,\  .        J        ft       ikzCOSa     /  _ri\ 

v«(P)  -  -A  sin  6  e      v(p,n-^). 


In  a  similar  way  we  conclude  that 


(63) 
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<i/-o\  .        „       a       ikZCOSO       /_    Ai\ 

w-(P;  -  A  cos  8  e       w(p,P) 


vAP)     -  -A  cos  &   e       w(p,n-9F), 


where  w(p,C^)  is  the  function  defined  by  (2U)  of  Part  II. 

Upon  substituting  (82)  and  (83)  into  (79)  we  obtain 

E^'(P)  «  A  a  sin  e  J  p'^cos  ^fv(p,(?)  -  v(p,Ti-^)| 
(3h)  *   (ap)-^(p2-a^)^/2g.^r(p^^)  ^  v(p,ni20j 

p"-^sin^  [w(p,0)  -  w(p,n-0)J 


-  Aacos9cosa 


-  (ap)"^(p2-a^)^/^co8  0  pp,0)  ♦  w(p,n-;^)l 

f'^'(P)  -  A  a  sin  e  ^  p'^sin  0|v(p,0)  -  v(p,n-0)| 

(85)  -  (ap)"^(p^-a^)^/^cos/)p(p,C()  *  v(p,n-0)~| 

"  COS0  w(p,0)  -  w(p,n-^) 
'in(^    w(p,0)  +  w(p,n-<jf) 


+  Aacosdcosa  <  P 


,   .-1,  2  2x1/2  . 
♦  (ap;  (p  -a  )   sin' 


,d' 


(86)     E  (P)  -  -A  COS  ©  sin  a  (  w(p,0)  +  w(p,7t-0)  )   . 


z 


The  expressions  for  the  components  of  the  electric  field  are  primed  to  indicate 
the  supresaion  of  a  common  e^      factor.  If  we  now  use  the  fact  that  the 
cylindrical  components  of  E  are  given  by 
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F,    -  E   COS  0  +  E   sin  0 

d»       d*     _,    d'     _» 
E^   -  -E^  sin  0  ♦  E   cos  (^ 

we  see  that  (8h),  (85)  and  (86)  a^ree  with  the  sum  of  (?5)  and  (26)  of  Part  II 
which  is  the  asjmptotic  expansion  of  the  exact  solution  of  the  vector  boundary 
value  problem. 

10.  The  field  near  caustics 

The  caustic  of  a  given  family  of  rays  is  defined  as  that  surface  or 
curve  which  is  everywhere  tangent  to  the  ray  system.  Thus,  the  surface  of  a 
diffracting  body  is  a  caustic  of  the  family  of  diffracted  rays  resulting  from 
a  given  field  being  incident  on  the  body.  In  addition  to  the  diffractin?  bodies 
in  the  problems  we  have  considered  an  axial  caustic  (the  polar  axis)  is  formed 
in  the  sphere  problems. 

In  accordance  with  optical  experience  and  as  we  see  from  equations  (11) 
and  (2U)  on  a  caustic  the  field  constructed  by  the  present  method  becomes  in- 
finite, A  finite  value  for  the  field  on  the  caustic  and  a  corracted  value  for 
the  field  near  the  caustic  can  be  obtained  by  a  slight  modification  of  our 
method. 

Let  us  first  correct  equation  (11)  to  hold  on  or  near  the  body,  lb  do 
this  we  make  use  of  a  canonical  problem,  the  cylinder.  Equation  (5)  of  Part  II 
was  obtained  from  (U)  by  using  the  Debye  expansion"-  -I  of  the  Hankel  functions 
H^  (kr)  and  H   (kp)  appearing.  However,  when  we  are  near  or  on  the  body 
T  -^  a   and  since  v  ^  ka  this  expansion  is  not  appropriate.  In  this  case  a  valid 
approximation  is  the  Hankel  expansion*-  J  of  the  Hankel  Function.  Thus,  we  can 
modify  ( 5)  to  hold  near  or  on  the  body  if  we  multiply  it  by  the  ratio  of  the 


1^ 
To  be  precise  (5)  represents  the  leading  terra  in  the  asymptotic  expansion  of 

the  field  as  long  as  r  /  a.  However,  when  r  is  near  a  this  leading  term  has  no 

value  for  practical  computational  purposes. 
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Hankel  expansion  of  H^     (kr)  to  the  Debye  expansion,     tquation  (8?)  below  gives 
th«»  ratio  R/r     of  the  Hankel  to  the  Pebye  expansion  for  r  near  a, 

^1/2,1/3  -1/2.1/6,  2     2a/U        f-w  2     2a/2  -1/    /  n     in 

(87)  /     "  2       6  '^Tt     '   k  Mr  -a  )  '   exp  -ik(r  -a  )  '     +  iv;  cos     'a/r)-  =^ 

k  \  IkCa-r)   *  q^(ka)^/3^  (kr) 


We  now  postulate  that  on  or  near  the  body  (11)  should  pe   modified  by 

multiplying  by  R^  ,  This  requires  that  the  quantities  r  and  a  appearing  in  the 

above  be  given  meaning  for  an  arbitrary  convex  body.  The  radius  a  we  interpret 

in  the  same  way  as  we  did  in  the  expressions  for  a  and  D  ,  To  interpret  r 

ID     m 

2  2  1/2  2   2  1/2 

notice  that  s  -  (r  -a  )  '   so  that  if  we  set  r  ■  (s  +  a  )    then  r  is  defined 

in  terms  of  geometrical  quantities. 

Notice  also  that  in  case  the  source  is  located  on  the  body  we  can  modify 
equation  (5)  so  that  it  remains  valid  for  p  «  a  by  multiplying  it  by  R/;   with 
r  replaced  by  p.  This  admits  of  an  obvious  generalization  to  a  convex  body  on 
which  a  point  source  is  located. 

The  derivation  of  the  asymptotic  expansion  of  the  field  on  an  axial 
caustic  has  been  previously  studied  ^^  J,  The  method  consists  of  writing  down 
an  exact  solution  of  the  reduced  wave  equation  which  possesses  an  axial  caustic 
and  then  using  the  exact  solution  to  evaluate  the  field  on  the  caustic. 

In  cylindrical  coordinates  (x,p,0)  the  function 

(88)         u  ■  A  e       J  (kp  cos  5)  cos  n0 

is  an  exact  solution  of  the  wave  equation  for  arbitrary  5  and  A.  If  we  then 
asymptotically  expand  the  integral- ordered  Bessel  function  in  (68)  we  can  re- 
cognize (88)  as  describing  a  family  of  rays  which  possesses  the  x  axis  as  a  caustic 
and  makes  an  angle  of  n/2  -  6  with  it.  By  then  comparing  (88)  with  its  asymptotic 
expansion  it  can  be  easily  shown  that  to  correct  a  given  optical  field  which 
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possesses,   i.e.,   becomes  infinite  on,   an  axial  caustic  we  must  multiply  the  ex- 
pression  for  the  field  by 

(89)  (l/2)(2nkpco86)-^'^  secrkpcos6  -  (n+  i)  i  1   J  (kp  cos  5)    . 


;  kpcos6  -  (n+  ^)   ^   J^( 


Equation  (89)  pivas  the  correction  factor  on  or  near  any  axial  caustic.  If  we 
are  on  the  caustic  (p  ■  0)  then  we  must  multiply  the  optical  field  by  (89)  and 
then  let  p  tend  to  zero. 

Let  us  now  apply  (89)  to  the  case  of  back  scattering  from  a  sphere.  In 
this  case  p  »  r  sin  6,  cos  5  -  a/r,  n  -  0,  and  (89)  becomes 

(90)  (2nka  sin  6)^^^ 

2  cos  (  J  ) 

Thus  for  ©  "  0,   equation  (11)   of  Part  II  becomes 

(91)  u(P)  --  u  (P)  +  lim       iinlSHiHl) n(P). 

S     e-^    2  cos  (J)     "^ 

As  an  example  of  the  numerical  usefulness  of  the  p^ometric  theory  we 

use  (91)  to  compute  the  back  scattering  from  a  rigid  sphere  upon  which  a  plane 

ikz 
acoustic  wave,  u.  «=  e   ,  is  incident  and  compare  the  res-alt  with  a  computation 

based  upon  sximminp  the  exact  series  solution.  For  the  geometrically  reflected 
wave  u  (P)  we  use  the  leading  tern  of  a  previously  calculated  asymptotic  ex- 
pansion"- ""J 

(92)  u  (P)  ^   -   a(2r)~-^  exp|ik(r-2a)l  . 

Equation  (92)  is  the  far  field  result,  i.e.,  for  r  tending  to  infinity. 

If  we  follow  the  method  of  Section  6  for  the  case  of  an  incident  plane 
wave  we  obtain  for  points  in  the  lit  region,  but  not  on  the  axis 

(93)  u^(P)  ~  w(r,e)  -  iw(r,-e). 
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In  (93) 

(9U)  w(r,e)  - 


=SSS|SP  r.:  ...,[=■•<-..]■' 


exp 


(ik-ajjj)a  -^  |^  ~  ©  ~  cos"  (a/r)^ 


The  values  of  D  and  a  in  (9U)  are  given  in  Table  I,  Cas«  II.  In  the  following 
in     n 

we  will  use  only  the  node  corresponding  to  in  ■  0.  When  we  substitute  (92)  and 
(95)  into  (91)  and  let  r  tend  to  infinity  we  obtain 


(95)     u(P)  -  -a( 


where 


2r)"-'-  exp[ik(r-2a)l  F(ka) 


(96)     F(ka)  -  1  -  2A26"^/^q;V2(q^)(ka)^/\xpr-nq^6"^/^(ka)^/^8in  ^1 

expF^  +  ika(n+2)  ♦  iiiq^6"^/-(ka)^/5co8  ^l   . 


By  an  analysis  of  the  exact  series  solution  to  the  associated  boundary 

value  problem  it  is  easy  to  see  that  |F(ka) |     is  the  asymptotic  expansion  of  the 
seiries 

00  _  J„(ka)  |2 


(97)  iE(ka)r     -    -^2-        E  (-1)"  (2n*l)  ^p 


(ka)        '   n-o 


h^^^   (ka) 


In  Figure  8  the  functions  |E(ka)l  and  |F(ka)|  are  graphed  for  2  <  ka  <  10, 
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Appendix  I 

In  ordar  to  apply  •quation  (11)  for  u,(P)  to  the  sohere  it  is  necessary 
to  calculate  the  radius  of  curvature  p.  This  quantity  was  defined  as  the  non- 
zero principal  radius  of  curvatiire  of  the  wave  front  on  the  body.  Now  the  radius 
of  curvature  of  a  plane  curve  is  e-^ual  to  the  distance  from  the  curve  along  the 
normal  to  the  point  where  neiphborinp  normals  intersect.  Furthermore,  the  rays 
are  normals  to  the  wavefronts.   Therefore,  the  radius  of  curvature  of  the  normal 
section  cut  out  of  the  wavefront  through  P.  or  ?2   by  the  plane  of  Figure  5  is 
zero.  Since  the  principal  directions  are  orthogonal, p^  is  the  radius  of  curvature 
of  the  normal  section  obtained  by  intersecting  the  wavefront  with  a  plane  n  per- 
pendicular to  the  plane  of  Figure  S» 

Let  us  first  calculate  p^(P^),  The  trace  of  the  wavefront  on  the  plane 
TT  is  simply  the  point  P..  However,  p^  was  defined  as  the  limiting  value  of  the 
principal  radius  curvature  of  the  normal  section  of  the  wavefront  through  a  point 
P-  on  the  diffracted  ray  a»  P-  tends  to  P.,  If  we  again  use  the  fact  that  the 
radius  curvature  is  the  distance  frcan  the  curve  along  the  nonnal  to  the  point 
vdiere  neighboring  normals  intersect,  it  is  evident  that  the  limiting  value  of  p- 
is  as  shown  in  Figure  5.  This  is  the  distance  along  a  tangent  to  the  sphere 
through  the  point  P-  from  P,  to  the  point  of  intersection  of  the  tangent  and 
the  polar  axis.  A  simple  geometric  calculation  then  yields  (23)  for  p^(P,), 


A  normal  section  of  a  surface  is  the  plane  curve  cut  out  of  the  surface  by  a 
plane  containing  the  normal  to  the  surface.  Sinc-e  the  rays  are  orthogonal  to 
the  wavefront  the  normal  sections  through  P,  or  P„  are  planes  cut  through  P, 
or  Pp  containing  the  diffracted  rays  through  these  points. 
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Part  IT.  Asymptotic  Expansions  of  Exact  Solutions  of  Diffraction  Problems 

1.   Introduction 

In  Part  I  we  solved  diffraction  problems  by  the  geometrical  theory  of 
diffraction.  Now,  instead,  we  consider  such  problems  as  boundary  value  problems, 
the  field  being  required  to  satisfy  one  or  more  partial  differential  equations. 
We  consider  special  cases  which  can  be  solved  by  the  method  of  separation  of 
variables.  The  solutions,  which  are  well-known  infinite  series  of  product  solutions, 
are  useful  only  for  small  values  of  ka.  By  utilizing  an  appropriate  transformation, 
due  to  G,  N.  Watson,  we  are  able  to  evaluate  these  solutions  asymptotically  for 
large  values  of  ka.  In  some  of  the  cases  which  we  treat  this  evaluation  has  been 
carried  out  before  for  special  boundary  conditions.  In  others,  such  as  that  of 

the  finitely  conducting  circular  cylinder  hit  obliq\iely  by  an  electromagnetic 
wave,  the  asymptotic  evaluation  has  apparently  not  been  done  before. 

In  Part  I  the  results  of  the  geometrical  theory  of  diffraction  have 
been  shown  to  apree  with  the  present  results. 


2,   Scalar  diffraction  by  a  circular  cylinder 

In  this  section  we  will  formulate  and  solve  the  boundary  value  problem  of 
diffraction  of  a  cylindrical  wave  by  a  circular  cylinder  of  radius  a.  Then  we 
will  expand  the  solution  asymptotically  for  large  values  of  ka.  We  follow  the 
method  of  W,  Franz"-  -■ ,  which  we  apply  to  various  cases  not  previously  treated. 

We  consider  a  scalar  field  u(r,9)  which  satisfies  the  equations 

(V^  ♦  k^)u  -  -r''^6(r-p)e(e)         r  >  a 
(V^  +  k^)u  -  0  r  <  a 

lim  r'-'^Cuj.  -  i  ku)  -  0  • 

3>->00 

A  description  of  the  Watson  transformation  with  particular  reference  to  the 
early  literature  is  given  in  the  introduction  of  [6J , 
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We  will  consider  four  different  kinds  of  boundary  conditions  at  r  -  a.  Thay 

correspond  to  different  physical  problems  and  are  listed  below.  In  the  first 

three  cases  the  second  equation  above,  which  applies  for  r  <  a,  is  to  be  omitted. 


Case 


Boundary 
Condition 


Physical  Itescription 
Acoustic  Electromagnetic 


u  -  0 


u  ■  pressure 
soft  cylinder 


perfectly  conducting  cylinder 


II 


u  ■  pressure 
hard  cylinder 


u  -  H 
z 

perfectly  conducting  cylinder 


ni 


-ikZu 


u  ■  pressure 
absorbing  cylinder 


u  ■  E    or  H 
z  z 

absorbing  cylinder 


IV 


u 


rh 


u'i^\ 

1..    iL 

u       continuous 
r 

bu  ■  pressiire     r  >  a 

^1       ^ 

b^u  "  pressure     r  <  a 

or 

bu  1         ■    fb-,  ul 
L    Jr«a        I-  1  -"r^a 

u  -  nk-^H, 

+                      •• 

b    _  density  of  medium 
ET       density  of  cylinder 

u  ■  mJc,   H 

BT  "  J.. 

k£n 


r  >  a 


r  <  a 


r  >  a 


r  <  a 


By  expanding  the  diffracted  field  in  cylindrical  wave  functions  and  imposing 
the  boundary  conditions  we  find 


a2  - 


+  <X) 


(1)  u  .|^   e 


ime 


m— 00 


'^.(^^<) 


01 


hI^W,)  ,  r  >  a 


Here  r^  and  r^  are  respectively  the  larger  and  smaller  of  the  quantities  r  and 
P. 

The  operator  i^  ,  which  is  determined  by  the  boimdary  conditions,  is 
defined  in  the  following  table 

riFix)     •  F(x) 


Case  I : 
Case  II : 


Af(x)  -  ^ 


Case  III:      /\F(x)  -  ^  F(x)  +  iZ  F(x) 

Case  IV:    S\F{x)   -  b  k^  /(k^a)F(x)  -  b^k  J^(k^a)  ^  F(x)  . 


In  Case  IV,  F(x)  denotes  a  cylinder  function  of  order  m. 

The  representation  of  the  field  as  given  by  equation  (1)  is  useful  only 
for  small  values  of  ka.  A  numerically  useful  result  for  large  values  of  ka  is 
obtained  by  employing  the  veil- known  Watson  transformation  which  leads  to 
an  asymptotic  expansion  of  the  field.  Letting  C  be  a  contour  encircling  the 
real  axis  as  shown  in  Figure   9  $   we  can  write  the  sum  in  (l)  as  the  following 
contour  integral 


(2)    u--^ 


/   iv(e-T,)  Hi^^(kr^) 


7 


sin  vn   (^H^^^(ka)L 


^^^[aH|,^^(ka)J^(kr^)-/lJ^(ka)H[,^^(kr^)]  dv 


When  V  is  replaced  by  -v  on  that  part  of  the   contour  C  for  which  Imv  <  0 


U2a  - 


V    plane 


Figurg  9;  The  path  of  integration  for  the  Watson  type  integral  repre- 
sentation of  the  field  diffracted  by  a  circular  cylinder. 


u    plane 


Figure  10 .•   The  paths  of  integration  for  the  integral  representation  of 
the  field  diffracted  by  a  sphere.  The  contour  C  encircles 
the  positive  real  axis,  while  P'  coincides  with  the  upper 
branch  of  C  for  Re  n  >  -  ^  and  is  as  shown  for  Re  n  <  -  i  . 
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(2)   becomep 


.(1) 


;p  ilH^     (ka)    L- 

Here  T  is  a  contour  in  the  upper  half  plane  running  parallel  to  the  real  axis. 
The   contour  t:!  can  be  closed  in  the  upper  half-plane  provided  the  point  r, 
e  lies  in  the  geometrical  shadow.     Then  the  integral  can  be  evaluated  by  com- 
putinp  the  residues  at  the  zeroes  v     of  -TLH^.     (ka).     This  yields 

co6v„(e-ti)     -^\  ^^^^       ,^^  (.) 

When  we  asymptotically  expand  the  Hankel  functions  appearing  in  (U),  this 
equation  takes  the  form 

Iexp[lk{(r^-.^)1''^.(P^V)V^}] 


u    ~     - 


E 


'2k(r^.a2)l/^(p^V)^/'^ 
exp[iv^©]   +  exp[iv^(2TT-6)] 


(^)  m  1  -  exp[2iiiv  "] 

•  exp  -iv^(cos~  (a/r)+cos     (a/p)^         -^ Pj-y 


AHj,^^(ka) 


L    "-  'J     |^r\H;-(ka) 


m 


Equation  (5)  contains  the  leading  terms  in  the  asymptotic  expansion  of  the 
field  in  the  shadow  region.  A  sliphtly  different  evaluation  of  the  integral  in 
(3)  leads  to  a  sindlar  result  for  the  diffracted  field  in  the  lit  region  plus 


Equation  (3)  is  exact  in  Cases  I,  II,   and  III.     In  Case  IV,   however,  where 
the  operation    f\.   is  a  function  of  v,  other  integrals  also  appear.     In  the 
case  of  opaque  obstacles   (e.?.,  lossy  or  highly  conducting  ones)  those  integrals 
are  negligible   compared  with  that  retained  in   (3).     Consequently  we  neglect 
those  integrals  in  our  further  developments. 


uu  - 


additional  terms  which  represent  the  geometrical  optics  field  u  .  This  evaluation 
is  carried  out  for  the  field  diffracted  by  a  sphere  in  Section  3.  The  result  (5) 
has  been  compared  with  the  geometrical  result  in  Part  I,  Section  h» 

The  values  of  v^  and  of  nHj^^^(ka)  r|^  QH^"'"^(ka)|'-'-  appearing  in  (5)  depend 

upon  the  boundary  conditions.  In  each  case  v  has  the  form  v  »  ka+q  (ka/6)  '  e  "'  , 
The  quantity  a  is  the  root  of  an  equation  given  in  Table  I,  page  Ik   .  The  function 
A(x)  appearing  in  this  table  is  the  Airy  function  defined  by 

00 

(6)        A(x)  -  I  coE(r^-xr)dr  . 

In  Case  IV  the  results  in  the  table  have  been  further  simplified  by  assuming 
|k-|  »  |kl.  The  factor  AH^^^^Cka)!  —  n^H^^^^Cka)!"  is  given  in  the  last 
column  of  the  table. 


3»    Scalar  diffraction  by  a  sphere 

We  will  now  treat  the  diffraction  of  a  spherical  wave  ly  a  sphere  of  radius 
a,  using  the  same  techniques  as  were  used  in  the  previous  section.  The  formulation 
of  the  problem  is  the  same  as  in  that  section,  the  delta  function  and  the  radiation 
condition  being  modified  to  correspond  to  the  present  three-dimensional  case.  The 
source  is  located  on  the  polar  axis  at  r  ■  p,  ©  =  0.  The  symmetry  of  the  problem 
ensures  that  the  solution  is  independent  of  the  azimuth  0» 

Upon  expanding  u(r,©)  in  spherical  wave  functions  and  applying  the  boundary 
conditions, we  obtain  the  solution 

(7)      „  .  g  I    (^n.l)P„(co*)ht"(.r,)[J„(^^,)-  ^j[]'!      -'"(kr,)]  . 

n  ' 

The  functions  h^'"  (x)   and  j^(x)  are  the  spherical  Hankel  and  Bessel  functions 


n+fr  ~"         n+: 


"n'W'/s    t>^-         V''-/h    \_1^-^ 
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The  operator    /\     is  the   same  operator  that  appeared  tn  the  cylinder  solution.     How- 
ever, in  Case  IV  the  Bessel  functions  in    ii     are  to  be  replaced  by  the  spherical 
Bessel  f'anctions  and  the  function  F(x)  is  a  spherical  cylinder  function. 

We  again  use  the  Watson  transformation  to  transform  (?)  into  a  form  which  is 
numerically  useful  for  larpe  ka.  Let  C  be  a  contour  encircling  the  positive  real 
axis  as  shown  in  Figure     10    ,     The   sum  in  (7)  can  be  written  as 

.        [       (2ti*l)Pj-cose)h^^^(kr^)  [l  (kr^)  Ah^^^ka)-  Ai  (ka)hf^^^kr^)] 

(8)  u--^  ti Si iJzli 1 ^^ ti ti ^dn' 

*^"    J^,  sin  nn  hj^^^ka) 

t 
Upon  replacing  p-  by  -jjl-I  on  that  part  of  the  contour  C  for  which  Imn  <  0  we  can  re- 

write  the  above  integral  as  an  integral  over  the  contour  D  shown  in  Figure  10  • 

The  integral  in  (8)  can  be  evaluated  in  terms  of  its  residues  at  the  zeros  a     of 

IB 

O    h       (ka)  provided  the  point  (r,©)  lies  in  the  shadow.     Then  the  integral  becomes 

(2n^+l)P     (-cose)  Ah^^^(ka) 

(9)  u.^r    '—^ h(^^(kr)h(^^(kp)  ^5 !!L_^ . 

TL.  sxn  ti^ii  ti^  (x^  |L.nh^^^(ka) 

If  we  set  V     ■  n     ♦  75-  than  the  v     entering  (9)   are  asymptotically  (for  ka  large) 
m         m       £  in 

equal  to  the  v     occurring  in  the  cylinder  solution.     We  also  have 

Ah^^^ka)                           HH^^^ka) 
m m 

l-Ah^^^ka)        ""         l-nH^^^(ka) 

Therefore, upon  asymptotically  expanding  the  Hankel  and  Legendre  functions  appearing 
in  (9), we  obtain 


Afain  this  is  exact  in  Cases  I,  II,  and  III  and  represents  an  approximation 
based  upon  the  opacity  of  the  obstacle  for  the  other  boundary  condition. 
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a 

,onkrpsin9)  '  (r  -a  )  '  (p  -a  )  ' 


exp[iv^(2n-e)-^]_.exp[iv^6.^] 


(10)  n         1  +  exp[]2niv  ^ 


nH^^^(ka) 


•  exp  -iv^jcos"  (a/r)+cos"  (a/pfj   -^—2 — j— 


-1 


(2)      r  9  "1" 

The  values  of  H.  H^  (ka)  ^Ah  (ka)    and  equations  for  a  ,  from  which 

V   •  ka  +  q  (ka/6)  '  ^exp{±n/3)   can  be  found,  are  given  in  Table  I,  page  li^  . 

The  above  result  for  the  field  in  the  shadow  does  not  apply  on  the  caustics 
r  -  a  and  9  "  n  because  the  asymptotic  expansions  used  for  the  Hankel  and  Legendre 
functions  do  not  apply  there, 

D-Ol 

We  will  now  apply  a  method  developed  by  W.  Franz  '-  -"  to  compute  the  field  in 
the  lit  region  by  means  of  the  Watson  transformation.  To  carry  out  this  evaluation 

we  write  the  Legendre  function  P  (-cos  9),  which  appears  in  the  integral  over  the 

t 

contour  D  ,  in  the  fora 

P  (-cos  9)     -     e^"^  P  (cos  9)   -  2i  sin  \in  Q^^^(9). 

f*  H*  r^ 

We  then  split  the  integral  into  the  sura  of  two  integrals,  one  containing  e  ^'   (cos  9) 

(2) 
and  the  other  containing  -2i  sin  |xn  Q   (9),  The  first  integral  can  be  evaluated 

by  residues  to  yield  the  diffracted  field.  The  second  integral  can  be  evaluated 

by  separating  out  the  incident  field  and  then  using  a  saddle  point  approximation 

to  yield  the  reflected  field.  Upon  applying  this  procedure,  the  field  in  the 

illuminated  region  is  found  to  be 


hi  - 


U  '^  u      ■♦• 


exp[ivje+2n)-  ^]  +  exp|iv^(2n-d)>  ^] 


E 


aHf.^^(ka) 


(11)  »        1  *  «^t2i"V„:| 

•  exp  -iv  ("cos"  (a/r)+cos"  (a/p)  T   -5 ^U-r 

L  "^  -U  ^OH^^Cka) 

m 

This  result  is  not  valid  on  the  caustics  r  ■  a  and  6*0  for  the  reason  given 
above.  The  diffracted  fields  in  (ID)  and  (11)  have  been  compared  with  the  geo- 
metrical results  in  Part  I,  Section  6, 

U»    Electromagnetic  diffraction  by  a  finitely  conducting  circular  cylinder  at 

oblique  incidence 

Consider  a  circular  cylinder  having  propagation  constant  11.  located  in  a 
mediion  of  propagation  constant  k.  Suppose  a  plane  electromagnetic  wave  of  angular 
frequenpy  co  is  obliquely  incident  upon  the  cylinder.  Let  us  introduce  a  rect- 
angular coordinate  system  in  such  a  way  that  the  generators  of  the  cylinder  lie 
along  the  z  axis  and  the  direction  of  propagation  of  the  wave  is  parallel  to  the 
(y,z)  plane.  If  a  is  the  angle  between  the  direction  of  incidence  and  the  z  axis 
then  the  direction  of  propagation  d,  is  given  by  dn  "  (0>  sin  a»  cos  a).  Now, 
if  the  wave  is  linearly  polarized  and  has  an  electric  vector  making  an  angle  'j'  "  ® 
with  the  (y,z)  plane,  it  can  be  written  as 


•  ^  i^(y  sin  a  +  z  cos  a) 
in  9  e   "' 


E   ■  A  si 

X 


ik(y  sin  a  +  z  cos  a) 
E   "A  cos  ©  cos  a  e   "^ 

y 


,     0  .     ik(y  sin  a  +  z  cos  a) 
E   ■  -A  cos  0  sm  a  e 
z 
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and 

,,             Ak  Q     ik(y  sin  a  +   z  cos  a) 

H       ■ cos  9  e       "^ 

X  [MO 


„  Ak        .      „  ilc(y  sin  a  +  z  cos  a) 

H       ■     —     sin  9  cos  a  e      ^"^ 


Ak     .      „      ,  ik(y   sin  a  +   z  cos  a) 

H       ■ sin  9  sin  a  e       "^ 

z  VM 


In  order  to  simplify  the  subsequent  analysis  we  break  up  the  field  into  two 
partial  fields  f'",  h''  and  E^,  H^  defined  by 


X         "X 


-y      -y  h:  '  -  H 


and 


z       z  £-„    "  H    -  H'  '  -  0 

A      y     z 


(2)  e(2)   .  ^(2)   .  ^(2) 


Because  of  the   symmetry  properties  of  Maxwell's  equations  it  is  sufficient 

to  consider  only  the  problem  for  which  the  incident  field  is  given  by  E   ,  H   , 

2     2 
This  is  due  to  the  fact  that  if  we  inter chan^  E  and  H  and  replace  ^  by  -k  /u  ji 

«• 
Maxwell's  equations  are  left  invariant   ^     Therefore, to  obtain  the  field  correspond- 

2       2 
ing  to  the  incident  field  E   ,  H     from  the  field  corresponding  to  the  incident  field 

1      A  2     2 

E   ,  IT"  all  we  need  do  is  interchange  E  and  H,  replace  ^i  by  -k  /w  n  and  then  replace 


,  .        AsinS      k 


■5 — - — ■  ~~" ' 

In  the  exterior  medium  we  take  Maxwell's  equations  in  the  form  7  x  E  •  ijiccH  and 

-ik^ 
y  X  H  -  — —  E.     In  the  interior  of  the  cylinder  |a  should  be  replaced  by  ^,   and 


k  by  k-. 
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To  calculate  the  diffracted  field  due  to  the  incident  field  E  ,  H  we  first 
introduce  polar  coordinates  (p,9^)  in  the  {x,y)   plane.  Next  we  define  two  potential 
functions,  n{p,^)   and  v(p,90,  both  satisfying  the  equations 


(72*k'2)   O   .  0 


p  >  a 


i^'*^')  o 


0  <  p  <  a. 


t  2         2       2     1/2 

Here  k'   «»  k  sin  a  and  k-.  »  (k^  -  k  cos  a)  '    ,     In  terms  of  these  potentials  we  de- 

fine  the  vectors  E  (p,^)  and  H  (p,0)  in  the  region  p  >  a  by 

r 


(1?) 


\    V 


ikcosg     3u 
r~    ^ 


k     u 


ino) 


8v 
5r 


/       t 

H 
r 


(13) 


V 


ik       Bu  .,  3v 

-ryr     +    ik    COSO       

liwr     3^  ar 


ik        $u  ik  ccsg       3v 

M.W      5r    *  r  a5? 


k       V      . 


Ihen  the  vectors  E  «  E'(p,gf)e^^^'^°^^  and  H  -  h'(p,0)8^'^'^"^°  are  solutions  of 


Maxwell's  equations L  J, 


Af  an  example  of  this  representation,  the  incident  field  E  ,  K  can  be  re- 
presented in  terms  of  the  potentials u.  and  v.  defined  by 


(DO 
(15) 


u.(p,0)   -  -(k')"'^A  COS&  sina  e^*''^ 


v^(p,^) 


50 


In  the  interior  mediviiti  the  expressions  for  E     and  E^^  are  obtained  by  replacing 

I  t 

M-  by  (JL,   iii  the  above  equations.     The  values  of  E     and  H     are  obtained  by  replacing 

It  t  t 

k     by  k-.     The  expressions  for  H     and  R^  are  obtained  by  replacing  n  by  |x.   and  k 

by  k.    In  the  terms  multiplying  the  derivatives  of  u  and  leaving  the  terms  multiply- 
ing the  derivatives  of  v  unchanged. 

When  E  and  H  are  represented  in  the  above  manner,  the  continuity  of  the  tan- 
gential components  of  E  and  H  across  the  cylinder  imposes  conditions  on  the  poten- 
tials u  and  V.     These  conditions  are 


k'^ul  »  k^^ul 

Jpaa*  -Ip-a- 


p*a- 


k'^vl  -  k'^vl 

Jpaa+  -J 

ikcosa     3u        .        3v 


p»a+ 


ikcosa     3u       ,  3v 


p=a- 


ik       Su  ikcosa     3v 

IKa)      3r  a  d5? 


^ik^     -> 
1     3u 


p-a+ 


(JLji-j)      ^ 


ikcosa     3v 
a  3^ 


p»a- 


We  now  write  u  -  u^   +  u     and  v  »  v,   +  v     where  both  u     and  v     satisfy  the 
is  is  s  s  ' 

radiation  condition.     Next  we  expand  the  field  in  cylindrical  wave  functions. 
Upon  taking  account  of  the  form  of  u.   and  v.    and  the  radiation  condition  we  ob- 


tain for  this  expansion 


00 


(16) 


AcosGsina    r— 
u -^ ^ 

k  n— 00 


00 


J^(k  p)  -  a, 


„H^^'(.'p)] 


in0 


Acosesina    r—        .       .  ,.  '    v  ^in0 
k  n"-oo 


p  >  a 


0  <  p  <  a 
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Acosesina  r—      u(1)m  '  \  in0 
k      n«-c» 
(17) 

AcosGsina  r—    ,   t  /i  '  ^  in0 
k      n»-  m 

The  constants  a  ,  b  ,  c  ,  and  d  are  to  be  obtained  by  applying  the  four  boundary 

conditions.  Upon  doing  this  we  find  that  the  constants  a  and  c  have  the  values 

n  n 


J  (k'a) 
n 

a     ■ 


'  .,.  »    .  (1)..    . 


"      AiiPo^lvi^    J„(kia)       k'a    H^l^(k'a)    / 

(IB)  r4  -^  ^'^  '  '"^'''^ 


k       H^k^a     J„(lc[a)       k'a     J^(k'a)/     V(k'a)^       (k^a)^- 


2         2 

n     cos  a 


(19)         c 


_     2  k  n  cosa  1  /      1         _  1 

""^M^n(k'a)''      ST^^OT^     V(k'a)^  "     (k^a)''^ 


where 


^^       J^(k^a)         ,       H^-(ka)\/k^       ^         Jn^y)  1       «n 


I ,  » 


(1)V.^^^    /.2  j'(,'    )  ^       Hl^^(ka) 


Hk^a     J^(k^a)       k'a    7^^J7I)~ J  [J^  V^^[^       J^i\^)       k  a     H^^\k'a) 

/      1  1      \^      2         2 

-  /       ,     g     -     —7—7         n     cos  a     . 

Uk  a)'^  (k^a)V 

[13I 
These  results  agree  with  those  given  by  J.R.  Wait*-     ■■  who  obtained  them  without 

making  use  of  potentials. 

Let  us  now  expand  these  results  asymptotically  for  large  values  of  ka.     First 

we  observe  that  the  term 
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tends  to  zero  as  ka  and  k,a  become  infinite*.     Therefore  c    -v^  0  and 

n 

a  '^  1^1  I —  -tr    — ^ — r  J —  ~  — r  1 —  ) 

k"  jJL^k^   '^n^*"!^^   ^    n  ^^  ^^ 

When  the  above  value  of  a  is  inserted  into  (16)  the  resulting  function  u(p,9^)  can 
be  shown  to  be  the  solution  of  the  following  scalar  problem 

(7^  *  k'^)u     -     0  p    >     a 

(7^  +  k!^)u     -     0  0  <  p  <  a 

u     ■     u.   +   u 
i         s 

p->co  \     '^  / 

5u 

^-^    continuous  at  p   ■  a 

dp 


bu  -     b, n 

Jp-a*  -Jp-a- 


Here  b  and  b-    are  given  by 


(21)  b     -     k'^n  /  k^ 


(22)  b^  "     y^\/k 


1 


If  thi5  term  is  retained  it  can  be  shown  that  it  will  not  influence  the 
leading  term  of  the  asymptotic  expansion.  Consequently,  we  set  it  equal  to 
zero  at  this  point  to  facilitate  the  subsequent  analysis. 
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Tt  is  to  be  noted  that  in  the  special  case  in  which  k.  becomes  infinite,  which 
occurs  when  the  conductivity  of  the  cylinder  becomes  infinite,  (20)  simplifies  to 

J^Ck'a) 

In   this  case  the  boundary  condition  on  u(p,0)  becomes  u(a,0)  "0, 

We  now  proceed  in  exactly  the  same  way  as  we  did  in  Section  2.  For  a  point 
p,  0  in  the  shadow,  we  find  that  for  large  values  of  ka  and  k^a,  provided  |k-  |  »  |k| 
u  is  asymptotically  given  by 

(23)         u  A^ -(k  )~  A  cos  6  sin  a  w(D,Cf)  +  w(p,n-0)  . 

Here  w  is  defined  by 

D^  exp[ik'(p2-a2)^/^  *  iv  (0  -  cos'^(|))] 


(2U)         w(p,0)  .    JZ   77-TTT7? 


m      (p  -a   ) 


1  -  exp[?T,iv    1 


The   values  of  D     and  v     can  be  determined  from  Table  I,  Case  IV  with  k  replaced 
mm  ' 

by  k   ,   k-   replaced  by  k.    and  b  and  b^    given  by   (21)   and   (22).     A  similar  result 
can  be  obtained  for  points  in  the  lit  region. 

We  now  calculate  the  vector  F     by  substitatinp   (23)  into   (12)   and  keeping  only 

I 

the  hif^hest  order  terms  in  k.  Upon  using  the  fact  that  v  -^  k  a,  we  obtain 


'    /  2  2vl/P  -1 

E  -^  A(p  -a  )  '"p  cos  9  cos  a 


p(p,0)  +  w(p,n-0)l 


(25)     Fgj  ~  Aap~  cos  &   cos  a  w(p,gf)  -  w(p,n-5f) 
F  ^  -A  cos  e  sin  a  w(p,0)  +  w(p,tt-0)   , 

The  electric  field  E  in  the  shadow  is  obtained  py  multiplying  (25)  by  e     °, 

2   2 
The  corresponding  result  for  the  incident  field  E  ,  H  can  be  obtained  from 

this  result  by  performing  in  (25)  the  substitutions  described  before.  The  result  is 


Sh  - 


E^    ^    A  a  p"-^   sin  9  v(p,^)  -   v(p,r.-0) 

E^   .-  -ACp^-a^)-*-/^  p"-'-  sinerv(p,(?)  +  v(p,i-0)j 


E      ^    0      . 
z 

The  function  v(p,0)  in  (26)   is  given  by  (2)4)  with  D^^  and  v^  determined  from  Table 
I,  Case   IV  with  b  and  b,    given  by 

(27)  b     -     k'^/n  b^     -     ^l^/^h.  ' 

These  values  cf  b  and  b^  are  obtained  by  performing  the  required  substitutions  for 

2 
H  and  y-.   in  (21)  and  (22),  Then  the  common  factor  -co  is  cancelled  from  both  of 

the  resulting  equations  since  only  the  ratio  bA^-i  occurs  in  the  problem.  The 

electric  field  E  in  the  shadow  is  now  given  by  multiplying  (26)  ty  e      , 

The  total  field  in  the  shadow  is,  apart  from  the  factor  e  ^^   ^°,  the  sura  of 

the  fields  (25)  and  (26)  due  to  the  incident  fields  E^,  H'^  and  E^,  H^,  The  field 

in  the  illuminated  region  consists  of  a  geometrical  optics  field  and  a  diffracted 

field  of  the  same  form  as  that  found  above.  The  above  result  for  the  diffracted 

field  has  been  compared  with  the  geometrical  result  in  Part  I,  Section  9, 


?,   Electromapnetic  diffraction  by  a  sphere 

We  will  now  consider  the  problem  of  diffraction  of  a  linearly  polarized  plane 
electrcmametic  wave  by  a  sphere  of  radius  a.  By  usine  techniques  similar  to  those 
used  in  the  previous  sections  we  will  solve  for  the  field  and  evaluate  it  asympto- 
tically for  large  values  of  ka. 

First  we  introduce  a  coordinate  system  with  the  origin  at  the  center  of  the 
sphere.  We  orient  the  axes  in  such  a  way  that  the  direction  of  propagation  of  the 
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incident  wave  is  the  positive  z  axis  and  the  direction  of  the  incident  E  vector 

is  the  X  axis.     The  incident  electromatTietic  field  can  then  be  written  as 

E     «E     «H     -H     -Oand 
y         z         X         y 

(28)  E       -     A  e^'''' 

X 

{29)  H       .     ^    e^^''     . 

Y        m 

Upon  introducing  a  spherical  coordinate  system  we  find  that  the  components  of 

[121 
the  diffracted  field  are  given  by^    -• 

(30)  E^  .  ii«22  ^    f      i"(2„*l)P„(cose)  (o„(»r)-b„h("(kr)) 

n«l  \  -^ 


Tf         iAcos0       9 
e  kr        ~" 


3       r      00    ^n 


^['&4fSSiV->(v-)-„Hi"<v.r.)] 


(31) 


.       „,     ^       00       .n 
Aco^ 


SI 


E  ^''-^      ^' 


°^    Is    r       '   P";V  P  (cose)  (i  (kr)-a  h^^)(k^)^ 
in6     W    (-^       nXn+T)       n  V^'^n^  n  n  J 


^4tS,^:{i?j!^v=-)0„(-)-vi^'<-')] 


(32) 


0  sin 

2     00      .n 


36     n«l  ^ 


Here  a     and  b     are  defined  by 
n  n 


^-  -ix»k,a  >-  -Jx« 


'^ 


ka 


(33)  « 


n 


>-  -x«k,a  *-  -' 


x»ka 
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(3U)         b__  . i 

L-  -ixok,  a  '—  -"x^k 


x«ka 


For  large  values  of  x,|^    xj  (x)     is  asymptotic  to  xj^^(x)  and  ^    xh^(x)l 
is  asymptotic  to  xh^    '   (x).     When  these  results  are  used  in  (33)   and  (3li)  for 
ka  and  k_a  large,  these  equations  become 

^xk^3jka;j;(k^a)  -  ^kj^(k^a)3;(ka) 
"         Mkj^h^-^^(ka)j^(k^a)  -  n^kj^(k^a)hj^-'^   (ka) 

k^  -  4  • 

—  k,i  (ka)i   (k,a)  -  —  kj  (k,a)j   (ka) 
M.      l^n^      ^n     1  m      ''n^   1     "^n^      ' 

(36)  b    -  ^ 


n        —  jT 

iL  k,h^^^(ka)j'(k,a)  -  J:  kj   (k,a)h^^^ '(ka) 
\i      In  n     1  \u      ''n     In 


By  the  same  technique  as  was  used  in  the  previous  section,  it  can  be  shown 
that  for  ka  and  k^a  large  and  the  point  r,  0.  in  the  shadow,  the  following  asympto- 
tic relations  hold 

-i  " 
(37)  21     i"(2n+l)P^(cose)  fjjkr)  -  b^h^^^kr)")  -.  v(r,e)+  e       ^  v(r,-e) 

n"l  V 


(3P)  Y_    i"  ^1^  P^(cose)  CVkr)-Vn^^^''''0    "  ^^^^^'^jj^^'®)*  «  ^  v(r,-e)J 

.  IT 

(^^^  £,   ^"  HfHTITPn(^°^)    (^j^(kr)-a^h^^^kr))-(ka)-^[^(r,e)  *e  ^  ^  u(r,-©)]  • 

Here  u(r,6)   and  v(r,9)   are  defined  by 
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*  -im  1  +  expQiv  n^ 

The  values  of  v     and  D     in  the  expression  for  u(r,9)   are  to  be  determined  by 
mm  "^ 

the  entries  in  Table  I  for  Case  IV  with  b  =  n,  and  b,    ■  m  •     The  values  of  v     and 

J.    J.  m 

C  in  the  expression  for  v(r,9)  are  to  be  determined  from  the  same  table  by  setting 

2  2 

b  -  k  /|i.  and  b.   ■  k./|j_..  These  facts  can  be  seen  by  comparing  (35)  and  (36)  vd.th 

the  expressions  for  the  f\   operator  in  Case  IV  given  on  page  xlj  • 

Now  to  calculate  the  electromagnetic  field  we  insert  (hO)  into  (37)-(39)  and 

use  the  results  in  (30)-(32).  The  resulting  expressions  for  the  electric  field 

components  in  the  shadow  are 

-i  " 
(Ul)      Ej.  ^   Aar'-''  cos  (2ffv(r,6)  -  e     v(r,-e)  | 

_.  n 
(h2)      Eq  -  A(r^-a^)-'-/^  r"^  cos  gfp(r,e)  +  e   ^  v(r,-e)l 

-i  " 
(U3)      E^  ^  -A  sin  9rpr,e)  +  e   ^  u(r,-e)l  - 


In  these  expressions  we  have  dropped  lower  order  terms  in  ka.  The  corresponding 
results  for  the  field  in  the  illuminated  region  are  obtained  as  in  Section  3.  The 
above  results  are  not  valid  on  the  caustics  or  on  the  shadow  boundary.  These  re- 
sults have  been  compared  with  the  geometrical  results  in  Part  I,  Section  8« 
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